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SUMMARY 


A program was conducted to determine if a state of the art micro- 
machined silicon solid state flow sensor could be used to replace the existing 
space shuttle orbiter flow sensors. The rather aggressive goal was to obtain a 
new sensor which would also be a multi-gas sensor and operate over a much 
wider flow range and with a higher degree of accuracy than the existing sensors. 
Two types of sensors were tested. The first type was a venturi throat design and 
the second was a bypass design. The accuracy of venturi design was found to 
be marginally acceptable. The bypass sensor was much better although it still 
did not fully reach the accuracy goal. Two main problems were identified which 
would require further work. 


Introduction 


The Gas Mass Flow Sensor Proof of Concept Testing for Space Shuttle Orbiter Flow 
Measurement program was conceived as a means of demonstrating that a commercially available 
state of the art micro-machined silicon solid state flow sensor would meet the stringent 
environmental requirements of the Orbiter. And also, hopefully, exceed the capabilities of the 
existing flow sensors by having a multi-gas capability and an improved flow range, improved 
accuracy, reduced power consumption and drift rate. Because of certain unforeseen problems, 
which will be discussed later, the program quickly became more development oriented rather than 
proof of concept oriented. As a result the time scale of the program was greatly extended and, 
although the program was intended to be fully funded by NASA, it became necessary for both the 
Honeywell Technology Center and Carleton Technologies Inc. to invest substantially in order to 
ensure that the results of the program approached the original goals as closely as possibly. Both 
Carleton and Honeywell were prepared to make this investment because of their continuing 
commitment to the concept of producing a space qualified flow sensor with currently unrivaled 
capabilities. 

The original program called for the manufacture and test of two through flow venturi 
design units (identified as Unit 1 and Unit 2) in order to estimate unit to unit variation. Fortunately 
this was done because a significant difference was found between the units (this will be discussed 
in the section on unit 2) which invalidated the premise that a multi-gas sensor could be produced. 
Also, substantial noise was found with this approach as well as great difficulty in reaching the 
accuracy requirements. For these reasons Honeywell recommended that a different approach be 
adopted. Honeywell had developed at their own expense a bypass type of flow sensor for another 
application and based on their estimation of its improved performance in the areas of noise, 
accuracy and unit to unit repeatability proposed that it be used as a replacement for the venturi 
type units. Honeywell offered to loan this bypass unit (identified as Unit 3) to Carleton and NASA 
for test purposes and as a result NASA agreed for the original plan to be modified to conduct only 
spot checks on unit 2 and introduce unit 3 into the program in place of unit 2. Unfortunately, 
because of a chip manufacturing problem (micro-cracks which are the subject of further 
discussion later) much time was lost in obtaining a properly functioning unit 3 and so the full range 
of tests was not conducted but rather a sufficient subset. 

This report shows an approximate time line of the testing which was conducted and 
mentions some of the minor problems encountered which are not addressed elsewhere. A 
discussion of the test results for each of the three units tested are presented in turn. The actual 
results plots are contained in three appendices (B,C & D respectively) to avoid cluttering the 
report. Appendix A contains a red-lined test plan and a red-lined test procedure and a paper which 
was prepared based upon the preliminary test data and presented at the 1996 ICES conference. 
The test plan shows a test setup schematic and details the testing methodology which was 
adopted in order to conduct a thorough and rational low cost test program. The paper shows 
cross sections for the two sensor types and gives a description of the principle of operation of 
each type thus this information is not repeated in the body of the report. Appendix E contains the 
Honeywell final technical report which was delivered to Carleton in March 1996. Finally, the 
conclusions section summarizes Carleton’s findings and recommendations for further means of 
advancing the technology to a point where it can be utilized to supplant the extant sensors. 


A Short Summary of the History of the Testing Conducted 


As the program was originally conceived it was believed that Honeywell would deliver a 
fully calibrated unit to Carleton for testing. However, Honeywell’s flow standard proved to be 
incapable of handling the wide range of flows to be covered, as a result, Carleton was forced to 
run all calibration tests before any significant testing could be conducted. These calibration tests 
consisted of six or seven flow sweeps covering ambient nitrogen, hot and cold nitrogen, oxygen 
and helium. Occasionally the cold nitrogen sweep would be dispensed with. After the calibration 
data was obtained a significant amount of manipulation was required to transform the data into a 
form easy for Honeywell to deal with. This calibration effort occupied between 4 and 5 days of 
work and as will be seen was repeated many times during the program. 

The original test program was scheduled to start in March 1 995 but was postponed to 
September 1995 because Honeywell experienced problems with the first batch flow sensor chips 
produced. These produced a large orientation error and Honeywell believed that by reducing the 
size of the chips by an order of magnitude the error would be significantly reduced. 

The first unit was delivered at the beginning of September and preliminary testing was 
started to try and get any bugs in the test setup worked out. The unit was observed to be very 
noisy. Also, the unit temperature sensor was observed to respond to unit body temperature rather 
than gas temperature. This latter fact caused an immediate change in the DOE structure as it had 
been planned to conduct tests in which the gas temperature and the unit temperature were 
different. Because it was essential for accuracy for the unit temperature sensor to give an 
accurate gas temperature reading it was necessary to only run tests in which the unit body and 
gas temperatures were identical. 

Calibrate unit 1 . 

Second unit delivered in middle of September. Unit observed to be extremely noisy. 

Calibrate unit 2. , 

Visit by Honeywell to review test setup. Agreed that unit 2 be returned to Honeywell to 
improve noise. 

Commenced testing on unit one. Observed that voltage out put not operating. Unit 
returned to Honeywell for rework. 

Unit 1 returned to Carleton with working voltage output. Commenced unit 1 calibration 
and observed that voltage output flipped to zero at maximum flow. Unit 1 returned to Honeywell 
for rework. 

Unit 1 returned to Carleton. Checked voltage output working correctly and noise improved 
by new software code. 

Honeywell by studying calibration data determine that there is a problem with unit 2 in that 
helium gas flow error performance is different from unit 1 by approximately a factor of 2. This is 
believed to be a function of the flow channel. 

Carleton flow measurement equipment damaged by improper use by new technician. 
Equipment returned to vendor for repair (November 95). Repaired unit returned end of December. 
Equipment recalibrated and accuracy checked. 

Unit 1 calibrated and the first meaningful test data begins to be collected. 

Honeywell state that flow channel problem with unit 2 is insurmountable and offer loan of 
bypass unit (unit 3) as a potential solution. 

Unit 1 DOE tests run (March 1996) 

Unit 1 flow, pressure, temperature, gas type sweeps run. 

Unit 1 step response tests run. 

Unit 3 received and calibrated. 

Preliminary testing shows pressure sensitivity and lack of repeatability. Unit returned to 
Honeywell. It is discovered that flow sensor chip is cracked (according to Honeywell caused by a 
flaw in manufacturing process). Inspection shows many other chips from same batch also 
cracked. Honeywell find one which is not and repair unit 3. 

Unit 2 calibrated. It is observed that the noise on this unit is substantially improved. A 
subset of flow, pressure, temperature and gas type sweeps is run as a spot check per agreement 
with NASA. The helium gas factor of two flow error is still present. 

Unit 3 returned to Carleton. 



Unit 3 calibrated. 

A subset of flow, pressure, temperature and gas type sweeps is run as a spot check 
taking into account the limited time left in the program. Some stand alone pressure sweeps are 
run at constant flow and a representative step response is run mid May 1996. Also, a number of 
temperature experiments run where the effect of small gas to unit body temperature differences is 
investigated. 

Finally, during the first two weeks of June some repeat flow sweeps are run to investigate 
the effect of zero drift/shift. 



Test Results Unit 1 


The DOE results data and plots are presented in appendix B. The figure and table 
numbers presented in this section refer only to appendix B. Calibration data is not included. The 
data is presented in approximately chronological order. 

Prior to commencing the main test program a preliminary flow test, using nitrogen at room 
temperature, was conducted to verify correct operation of the sensor and the test equipment. The 
results are shown in Figures 1 and 2. Figure 1 is Honeywell mass flow versus actual measured 
mass flow and Figure 2 is the error (as a percentage of reading) plot for figure 1 . These show that 
the sensor gives an output with a linear trend over the required flow range but with a maximum 
deviation of 19% of reading. It was determined that, apart from a narrow flow range, the deviation 
was significantly constant and thus could be attributed to a potentially correctable calibration error. 
Therefore the test program was continued knowing that absolute errors would be unacceptable 
but in the knowledge that average error could be improved by more accurate calibration. 

The first structured experiment (L8), which was designed to verify that humidity, ambient 
air pressure and X-axis orientation did not influence performance, was successful in that it did 
demonstrate that these factors had a negligible effect on the sensor output as expected. The 
results are shown in Table 1 (the actual DOE matrix is in figure 3 of the test procedure). The 
ANalysis Of VAriance (ANOVA) of the flow signal is in Table 2 and showed that 98.43% of the 
variation in flow was due to the flow changes within the experiment which indicates, not 
unexpectedly, that as a flow sensor the unit is very good. And also, that the other factors within 
the experiment, humidity, gas pressure, gas temperature, ambient pressure and X-axis orientation 
produced output variation of only 1.42% and so had almost no effect on the flow signal. The 
experimental error was only 0.15% of the variation within the experiment showing that the 
conditions were accurately set and that there were no significant effects from outside sources. An 
ANOVA was also performed on the flow error (as a percentage of reading) and is shown in Table 
3 and indicates that 31.8% of the variation in error was due to flow, 6.4% due to gas pressure, 
17.1% due to gas temperature, 6.0% due to ambient air pressure, 25.5% due to humidity and 
3.8% due to X-axis orientation with 9.4% experimental error. None of these effects were 
statistically significant and it was hypothesized that the apparent effects, particularly the large and 
unexpected effect of humidity on the error were due to the use of a saturated experiment and to 
the difficulty in maintaining consistent gas and unit temperatures. To verify this a full factorial L4 
experiment was run with flow and humidity as the only factors. The results are shown in Table 4 
and the DOE matrix in Table 5. The ANOVA of the flow signal for this DOE is presented in Table 6 
and shows that the flow changes accounted for 99.9% of the variation in flow and humidity for only 
0.05%. The main effects plot is shown in Figure 3 and shows that changing relative humidity from 
25 to 85% has an imperceptible effect on flow rate. The ANOVA for the flow error is given in Table 
7 and shows that error variation was 49.5% due to flow, 10.7% due to humidity which again is not 
statistically significant and 39.9% due to experimental error and other factors not included within 
the experiment. 

The results for the L18 experiment are shown in Table 8. Because the flow levels 
selected for nitrogen and oxygen were beyond the capability of the sensor for helium, the helium 
runs were conducted at flow levels one-tenth of the nitrogen levels. As a result, to analyze the 
data, the helium results had to be coded (multiplied by 10) to make them comparable with the 
nitrogen and oxygen. The flow signal ANOVA (Table 9) showed that 99.43% of the flow variation 
was due to flow changes within the experiment with no statistically significant effects from other 
sources. The flow error ANOVA (Table 10) showed no statistically significant effects, and a 
residual of 55.6% which indicates that there was more variation in flow error due to experimental 
error and outside factors than there was due to the factors within the experiment. Of the five 
factors within the experiment, the gas pressure and the flow level, were the main contributors to 
flow error. The main effects plots (Figure 4) show that the sensor output is significantly linear with 
flow and virtually not affected by the other factors. Figure 5 shows the interaction plots which 
seem to indicate the presence of many interactions. 

The LI 8 experiment was found to have been poorly designed because the flow levels 
selected were not equally spaced and because the interactions are uniformly distributed across all 
columns which gave the indication of many interactions. The lack of equal flow level spacing was 



overcome by analysis of the logarithmic flow value. A series of five L4 experiments were 
subsequently necessary in order to establish which interactions were real and which were 
spurious. The interactions indicated by the LI 8 were between gas pressure and gas temperature, 
between gas temperature and Y attitude, between gas pressure and Y attitude, between gas 
pressure and gas type and between gas temperature and gas type. The results and test matrices 
for the five L4 experiments are shown in Tables 11 through 15. The ANOVA results for the five 
DOE’s are shown in Tables 16 through 20. The main effects plots are shown in Figures 6, 8, 10, 
12 and 14 and the interaction plots in Figures 7, 9, 11, 13 and 15. It was recognized that DOE 3D, 
the gas type and gas temperature interaction DOE, had been incorrectly run by using helium flow 
rates which were out of range so it was re-run and the revised matrix and results are shown in 
Table 21, the ANOVA in Table 22 the main effects plot in Figure 16 and the interaction plot in 
Figure 1 7. The interaction plots for the 5 DOE’s show that the only real interactions are between 
gas pressure and gas temperature and between gas pressure and Y attitude. Both of these 
interactions are relatively weak being about 1.4% and 4.5% of flow reading respectively but 
together make it almost impossible to achieve the target 3% of reading requirement. 

The gas pressure and Y attitude interaction is caused by the well known gravitational 
effect and was already identified by previous testing by Honeywell. The L4 results show that at low 
pressure (100 psi) the gravitational effect is almost zero but at high pressure (900 psi) is quite 
large such that when gas is flowing up through the unit the measured flow rate is approximately 
4.5% higher than when nitrogen gas is flowing down through the unit. 

The gas pressure and gas temperature interaction was previously unknown and shows 
what appears to be a gas density related effect such that at low temperature and high pressure 
(when gas density is a maximum) the unit measured flow rate output is about 1 .4% higher than 
when the pressure is low or the temperature is high. As will be seen in the section on Unit 3 test 
results this was proven to be not a real unit interaction but a spurious one caused by lack of chip 
thermal isolation combined with the inability to keep the actual gas temperature and unit 
temperature exactly equal throughout the test. The DOE result was real but the cause was 
unfortunately not understood at the time. The program should have been halted to really 
understand this effect and then, in the long run, much time would have been saved. 

Table 23 is a summary of the data found in the following figures 18 through 145. The 
Table shows clearly that average error and standard deviation is relatively constant throughout 
gas type and pressure and temperature ranges covered. This indicates the potential to improve 
the performance of unit 1 by calibration and filtering in the flow range covered by the Gdu signal. It 
was this information which gave confidence that the unit 3 results would be much better than unit 
1 and the realization that attempting to cover the flow range by merging different flow sensor 
methodologies was not the best way to proceed. 

Figures 18 through 145 are plots of the smart sensor flow signal output versus the actual 
measured flow for each of the seven laminar flow elements (LFE’s) necessary to cover the flow 
range. These graphs show the linearity and hysteresis performance for one particular gas type, 
pressure and temperature combination and are representative of the other data collected. 
Included for each of the 12 cases shown in Table 23 is a figure which shows a concatenation of 
the seven individual tests and a figure which shows the error as a percentage of reading. These 
graphs show clearly the problems associated with measuring flow rate over such a wide range 
and problems with this particular unit. 

Noise at the low and high flow ends of the flow range is seen. The low end noise is 
caused by poor signal to noise ratios and electrical noise on very low signals. The high end noise 
is caused by flow turbulence and the inability to create a single flow passage capable of 
maintaining laminar flow over such a wide flow and pressure range. 

On the low flow plots (below 0.03 lb./hr for nitrogen) with LFE 1 a phenomenon is 
observed in which a vertical line of data points lies off to the side of the main response line. This is 
an error with the Flo-Dynamics computer used to convert the LFE signals into actual mass flows. 
It is caused by noise on the delta p transducer at very low flows causing the transducer signal to 
become negative signaling reverse flow through the transducer. The computer, instead of clipping 
or zero limiting these spurious values, took the absolute value and thus reflected them back so 
they appear to be valid data in the actual flow range. 


The concatenated error plot shows that the concept of using two different flow 
measurement methodologies to cover the high and low flow ranges with a merging of the two 
methods in the middle is not capable of producing low errors in the merged area. This may be 
capable of being corrected by use of multiple lookup tables instead of curve fitted equations but a 
superior approach appears to be that which is embodied in the by-pass unit which uses only a 
single measurement methodology and avoids the merging of two separate range signals 
completely. 

Figures 146 through 159 show step response tests for the sensor for the particular gas 
type, gas temperature and pressure conditions covered by the flow sweep tests shown in Table 
23. The flow rate selected for the test was 8 Ib./hr for oxygen and nitrogen and 1 Ib./hr for helium. 
The sensor does not exhibit a typical linear system exponential time response, instead having a 
pure time delay of about 2.5 seconds before rising (or falling) instantly to the final value. The 
sensor uses a digital sampling and averaging technique to minimize the effect of flow noise and 
calculate the output mass flow signal, and thus changes the output once every 2 to 3 seconds. 
This low update rate is related to the use of non-optimized electronics in the unit and would be 
expected to be increased in a production unit by the use of a more powerful micro-processor. 

Because of the averaging method used to reduce the noise levels the time delay 
associated with changing signals makes it very difficult to distinguish between true hysteresis and 
time delay effects. The faster a sweep is run the wider the apparent hysteresis band due to the 
time delay. The best estimate of hysteresis is obtained by assuming that the narrowest band is 
representative of the hysteresis and is least affected by the time delay. Figures 160 through 163 
show a selection of flow sweeps in which the direction of increasing and decreasing flow has been 
shown. Figure 164 shows a segment of one of these flow curves with best fit straight lines for 
increasing and decreasing flow. The hysteresis is defined as half the difference between the two 
best fits and is approximately 0.7% of reading. The flow standard used has hysteresis of 0.1% of 
reading thus, the unit 1 hysteresis is estimated to be 0.6% of reading. 

The power consumption of the unit was measured by measuring the current on each of 
the three power inputs in turn. The results are tabulated below. 


Voltage level 

Current 

Wattage 

-15 V 

69 mA 

1.035 W 

+15 V 

64 mA 

960 mW 

+5V 

70 mA 

350 mW 


The total power consumption is about 2.35 watts. However, as the electronics are 
commercial and non-optimized it is expected that this figure will be capable of reduction in a flight 
unit. Also, several extra circuits related to heater power measurement and control are included in 
this design which would not be required in a bypass type of unit. 

The pressure differential across the unit was measured at various flow rates. The values 
obtained were quite large, much larger than the target maximum values. For example, the 
pressure drop for nitrogen at 150 Ib./hr is required not to exceed 28 psid and the unit measured 
value was 42 psid. It was theorized that the these large values were due to the filters and not to 
any inherent problem with the unit. This was tested by utilizing the existing Shuttle engineering unit 
flow sensor. It was equipped with the filters taken from unit 1 and the pressure differential was 
measured. The filters were then removed and the test repeated. The results are tabulated below. 




Drop without filter (psid) 

15 

1.0 

0.2 

75 

12.7 

1.6 

150 

49.2 

5.0 


Thus it was concluded that the bulk of the measured pressure drop was due to the filters 
and not the unit. Because it is not feasible to test unit 1 without the filters in place actual pressure 
drops cannot be measured until the filter problem is corrected but based on the current 


















observations it is concluded that the actual unit 1 pressure drop will be slightly less than the 
existing Shuttle flow sensor. 


Test Results on Unit 2 


The plots of the test data are presented in appendix C. The figure numbers 
presented in this section refer only to this appendix. Calibration data is not included. 

The data is presented in approximately chronological order. 

As already discussed in the introduction, due to changes in the test plan and limited 
testing time available only a few spot checks were run on unit 2. Just sufficient to see that the unit 
is functional and to obtain a rough measure of the performance. Figure 1 is a plot of unit 2 mass 
flow versus actual measured mass flow at a number of discrete points across the flow range for 
nitrogen gas at 200 psig and room temperature. As can be seen the unit is, like unit 1 , 
approximately linear from 0.1 to 8 Ib./hrand then deviates substantially. The deviation occurs at 
the point where the unit begins to transition from a differential flow measurement methodology to 
a heater power methodology. 

Figure 2 is an error plot for the flow sweep shown in figure 1 . Error is presented as a % of 
reading and the different values obtained from each of the seven laminar flow elements LFE’s 
used to cover the range are separately identified. As can be seen the unit is not very accurate with 
the range for reasonable accuracy being only 0.1 to 2 Ib./hr. 

Figure 3 is a flow sweep for nitrogen at 200 psig and 130 degrees Fahrenheit and is 
similar to figure 1 . The error plot for figure 3 shown in figure 4 indicates about the same degree of 
accuracy. 

Figure 5 is again similar to figure 1 but is for helium at 200 psig and room temperature. 
Notice that the unit 2 measured flows are shifted relative to the actual measured flows by about a 
factor of 1.53. This is the reason the type 1 sensor cannot be used as a multi-gas sensor. 
Although unit 1 was acceptable as a multi-gas sensor unit 2 cannot be calibrated for both nitrogen 
and helium simultaneously because of this factor of 1 .53 error. For sensors operating correctly, 
such as unit 1 , a series of sensor signals versus flow can be obtained as a function of 
temperature and gas type. At constant temperature the differences in these sensor flow 
signatures for different gases amounts to a simple zero shift which allows all of the signatures to 
be collapsed into a single calibration curve. For the unit 2 sensor this was not possible as there 
appears to be an approximately constant multiplication factor relating the different gas flow 
signatures. Figure 6 is the related error plot and highlights the problem by showing the error plot is 
centered about 35% (1-1/1 .53) instead of 0 as with the nitrogen error plot. 

Figure 7 is a flow sweep for helium at 200 psig and 40 degrees F. It shows the same 
trend as the helium at room temperature. Figure 8 is the associated error plot which is again 
centered around 35% instead of zero. 

Figure 9 is a flow sweep for oxygen at 900 psig and 40 degrees F. It shows a similar trend 
to the nitrogen plots with a fair degree of linearity up to the sensor transition point. The equivalent 
error plot is centered about zero as expected and shows reasonable accuracy from 0.05 to 9 
lb./hr. 


Test Results Unit 3 


ri 
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The plots of the test data are presented in appendix D. The figure numbers 
presented in this section refer only to this appendix. Calibration data is not included. The data is 
presented in approximately chronological order. 

Due to the short time scale remaining in the test program only a limited amount of testing 
was accomplished, but, it is hoped sufficient to enable an appreciation of both the problems 
associated with and the capabilities of the unit. No designed experiments (DOE’s) were run on 
unit 3, however, as a full range of DOE’s was run on unit 1 and as unit 1 and 3 used identical flow 
sensing chips, simply arranged in a different geometry, it is believed that no loss of information 
occurred. The environmental factors covered by the DOE's should have generated the same 
results when applied to the chips in unit 3 as were obtained with the chips in unit 1 . Unfortunately 
only one version of unit 3 was available and so no estimate of unit to unit variability was capable of 
being obtained. This is the major shortfall in the testing conducted on the bypass unit. 

One reason for the short time scale is that a lot of time was lost due to the fact that the 
first time unit 3 was delivered to Carleton by Honeywell it contained a cracked flow chip which 
apparently caused many anomalies related to non-repeatable results. In general, data collected 
with the cracked chip is not presented. The cracks are the product of a defective manufacturing 
process utilized in the manufacture of prototype chips and not representative of production units. 

The tests which were conducted on unit 3 consisted of a series of pressure, temperature 
and flow range sweeps. The sweeps were not continuous as were the sweeps which were 
performed on unit 1 but were a series of discrete data points collected over the flow range. The 
flow sweeps were run for different gases, (nitrogen, oxygen and helium) at different pressures and 
temperatures. A similar spectrum was covered to that which has already been discussed in the 
test results for unit 1 . 

One of the first points to mention, with regard to unit 3, is that it was not optimized for the 
flow range requirements as units 1 and 2 were. Its main flow channel diameter is 25 mm whereas 
in units 1 and 2 the throat is 10 mm. Thus, theoretically unit 3 should be able to handle much 
higher flow rates than the venturi type units. However, Figure 1 , which is a plot of the raw sensor 
signal Gdu versus flow, shows a characteristic curl over at a flow rate of about 170 Ib./hr. This 
represents the absolute limiting flow rate of the unit because beyond this point the relationship 
between flow rate and Gdu is multi-valued and hence indeterminate. So it appears that although 
theoretically this flow channel should handle flows up to nearly 1000 Ib./hr it is only capable of 170. 
This is an unexplained anomaly which requires further investigation. Although these particular 
results were obtained with the cracked flow chip the replacement chip also showed a very similar 
trend in this area with an apparent curl over at about 190 Ib./hr. 

The second most important point to mention is that the unit is highly sensitive to small 
differences in temperature between the gas and the chip on board temperature sensor. 
Unfortunately this effect was not recognized until very late in the testing and so all calibration data 
was collected with only nominally similar gas and unit temperatures. This will probably explain 
much of the difficulty in obtaining a good unit calibration. This effect was misdiagnosed originally 
as either drift, or pressure sensitivity. Ultimately it was realized that the apparent lack of 
repeatability or pressure sensitivity was totally due to these temperature differences which were 
not being controlled as tightly as necessary because the result of very small temperature 
differences was not understood. 

Figure 2 is a flow sweep for nitrogen at 200 psig and 100 degrees F. Figure 3 is the 
associated error plot and Figure 4 is an enlargement of the error in the range ±10%. As can be 
seen the unit is not calibrated properly because there is a distinct trend with error moving from 
negative to positive as flow increases, in other words a slope error. Also, because of the curvature 
at the low flow end it is obvious that the unit has an incorrect value of Gdu-zero. The range of 
flows where the unit remains in the approximately ±3% range is from 2 Ib./hr to about 80 Ib./hr. 
Figure 5 is a flow sweep for oxygen at 900 psig and room temperature and Figure 6 is the related 
error plot again showing good accuracy between 2 and 80 to 100 lb./hr. Figure 7 is a flow sweep 
for nitrogen at 200 psig and room temperature and Figure 8 is again the error plot showing good 
accuracy between about 1 and 150 lb./hr. Figure 9 is a flow sweep for helium at 200 psig and 
room temperature and Figure 10 is the error plot showing good accuracy between about 4 and 7 
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Ib./hr. Figure 11 is a flow sweep for oxygen at 900 psig and 150 degrees F and Figure 12 is the 
related error plot. Figure 12 is the same error plot with an expanded scale and shows good 
accuracy between 1 and 40 Ib./hr. Figure 13 is a flow sweep for helium at 200 psig and 30 
degrees F and Figure 14 is the related error plot and shows good accuracy between about 9 and 
22 lb./hr. Figure 15 is a flow sweep for nitrogen at 200 psig and 150 degrees F and Figure 16 is 
the related error plot. Figure 1 7 is the same error plot with an expanded scale and shows good 
accuracy between about 3 and 30 lb./hr. This data is summarized in Table 1 below. 

Table 1 - Flow Accuracy Summary 


1 Gas 

Pressure 


Accurate Flow Range 

IHZfi&StEl&HS 

200 

Ambient 

1 to 150 


200 

100 

2 to 80 


200 

150 

3 to 30 


900 

Ambient 

2 to 90 


900 

150 

1 to 40 

Helium 

200 

Ambient 

4 to 7 

Helium 

200 

30 

9 to 22 


What these flow sweep plots and the table show is that the unit is capable of covering the 
flow range over a range of pressures and temperatures for all three gases tested. All of the plots 
show the similar trend of incorrect calibration at the low flow end. This is caused by a value of 
Gdu-zero (the basic sensor zero offset) which is incorrect. The fact that this bias value was set 
incorrectly is a combination of two factors, firstly, as already mentioned, the unit is very 
temperature sensitive and the gas to unit temperature difference was not sufficiently tightly 
controlled during the calibration runs. Secondly, the unit was suffering from a burn-in transition : 
drift which was not recognized at the time as will be explained later. As a result the unit was not in 
a stable situation for calibration or for subsequent testing. 

Figures 17 and 18 show a representative step response for the unit. The demonstrated 
characteristic is, as expected, identical to that which was obtained from unit 1 and reference to the 
section on unit 1 test results will provide expanded comment. 

Figure 19 is one of the many pressure sweep tests which were initially run and which 
caused considerable worry because they seemed to show that the unit mass flow output was a 
function of pressure and/or also varied widely with time { the ±3% of reading band is also shown 
for reference). However, Figure 20 shows that while the gas and unit temperatures were nominally 
kept constant throughout the test the temperature difference actually changed with pressure 
(time) and an almost perfect correlation is seen between the apparent flow variation and the 
extremely small temperature difference variation. 

Figures 21 through 27 show the results of a series of tests designed to investigate this 
temperature sensitivity. For these tests the actual mass flow, gas pressure and unit temperature 
were rigorously held constant. Initially the gas temperature was set as closely as possible to the 
unit temperature. After a suitable stabilization period the gas temperature was given a slight 
positive disturbance followed by a negative disturbance and finally returned to the unit 
temperature value. The curves show large flow reading changes (up to 40% for a 7 degree F 
temperature change). The flow variations show good correlation with the temperature difference 
variations as shown in Table 2. These results show excellent agreement with the results of the 

Table 2 • Flow to Temperature Difference Correlation 


Pressure 

Flow Rate 

Correlation Coefficient 

200 

2 

0.6852 

900 

2 

0.9783 

900 

10 

0.9758 






































DOE which was run of unit 1 . The DOE highlighted an unknown pressure and temperature 
interaction and indicated that flow rate was more sensitive to temperature at high pressure than at 
low pressure. As can be seen here, at low pressure, the correlation is not as good as at high 
pressure. 

It was decided to determine whether this effect could be corrected by temperature 
compensation. Because the unit temperature sensor was mainly responding to unit body 
temperature changes instead of gas temperature changes the unit effectively had no 
compensation. In these tests, where unit sensor temperature was deliberately kept completely 
constant, it was definitely completely ineffective. Thus the 900 psig and 10 Ib./hr case was 
manipulated on a spread sheet and a temperature compensation algorithm based on the actual 
gas temperature was applied. Figures 28 and 29 show the before and after results for mass flow 
and flow rate error as a function of time. It can be seen that application of the algorithm was 
capable of reducing the 14% flow deviation to about ± 2%. Unfortunately, the coefficients used in 
the algorithm were much larger than expected and would probably be inappropriate for larger 
temperature variations. Thus the whole question of temperature compensation for the unit is 
something which does require further investigation. Honeywell have provided the curve shown in 
Figure 30 which shows that the temperature compensation algorithm with normal coefficients 
works quite well over a very wide temperature range for their standard production unit. 

Normally, when Honeywell manufacture flow units the chips are subjected to a proprietary 
hardening or burn-in process. However, because the chips used in units 1 through 3 were new 
prototype chips not manufactured by the standard process they were not subjected to this 
accelerated burn-in cycle. (These prototype chips were used because they were designed to 
reduce or eliminate the orientation error.) Honeywell’s experience suggested that natural burn-in 
could be accomplished by normal operational use over a number of days. However, it appears 
that in the case of these new chips that was not true and as the testing cycled the units through 
higher and higher temperatures the zero offset continued to change. This effect is illustrated by 
the final two figures 31 and 32. The initial flow sweep was run after calibration at a maximum 
temperature of 130 degrees F and showed a large negative error at low flow. The unit was 
subsequently tested at 150 degrees and then the ambient test was repeated on 6/7/96 and 
showed a large positive error (the zero offset of the unit was noticed to have changed by about 3 
milli-volts between runs). One week later (the unit was kept continually powered) the test was 
repeated again on 6/14/96 with almost no change indicating that the unit had stabilized and could 
now be accurately calibrated at low flow. The high flow calibration accuracy would still be 
dependent on obtaining proper temperature compensation. 


Conclusions 


New programs seldom move to completion as originally planned and this one was no 
different. Hindsight indicates that more preparatory work/testing should have been carried out 
before launching such an ambitious program. None-the -less, much work was accomplished and 
much information was gained so that although the program was only partially successful in 
reaching the original goals, enough was learned to enable some cogent observations and 
recommendations to be made. 

The Accomplishments: 

It has been determined that both the venturi and bypass types of sensor are capable of 
operating over the required gas type, pressure, temperature, orientation attitude and flow ranges 
required. Also it has been demonstrated that the units utilize little power, have a fairly low pressure 
drop and a relatively fast response time. It has been shown that the units are insensitive to relative 
humidity and ambient pressure. In terms of flow accuracy and noise the venturi type is marginal 
but may have the potential to be improved to a point where it could meet the accuracy goal for a 
specific gas type. The bypass type can definitely meet the accuracy target over very wide flow 
ranges and for various gas types. The unit’s orientation sensitivity in a gravitational environment 
has been established and measured. The units have been shown to be both robust, having 
survived many shipments between Carleton and Honeywell, and capable of operating over 
extended periods without maintenance. The importance and required avoidance of certain 
manufacturing defects, such as micro-cracks, has been highlighted. 

The Remaining Work: 

What has not been achieved is that no measure of unit to unit variation has been obtained 
for the bypass unit and as was seen with the venturi sensor type this could be crucial. Satisfactory 
temperature compensation has not been demonstrated for either sensor type. Although Honeywell 
have supplied data showing good temperature compensation for production units the data do not 
address either high pressure or very small temperature differences, which effects have been seen 
to be large. Further, the ability to satisfactorily thermally isolate the chip from the unit body, to 
allow proper temperature compensation, has not been demonstrated. Also, related to this feature, 
the effect of thermal gradients caused by gas to unit body temperature differences has not been 
investigated. Finally, it has not been demonstrated that the bypass type of sensor can 
accommodate the required flow ranges when the flow channel is sized according to recognized 
theoretical guidelines. 

Recommendations: 

Carleton and Honeywell are still committed to the long term goal of producing a space 
qualified flow sensor with capabilities which are not available from other products. However, being 
both pragmatic and averse to taking unnecessary risks Carleton believes that before launching 
into a full flight development and qualification program a further small development/demonstration 
program should be completed to address all of the items listed above in the remaining work 
section. Carleton therefore proposes that Carleton and Honeywell be funded by NASA to modify 
the two units (Unit 1 and Unit 2 which are NASA owned products of the current program) for 
further development. This modification would involve re-using the existing electronics but 
replacing the flow channels with properly sized bypass type flow channels fitted with actual 
production chips (not prototypes). The programmed modification would be designed to allow 
initially, investigation of chip thermal isolation, and subsequently, the use of a separated 
temperature sensor located in the flow path as a fall back position, should sufficient thermal 
isolation for satisfactory temperature compensation prove impossible. 

In addition, to further lower risk, funds should be made available for the design, building 
and stand alone testing of a set of prototype flight type electronics. 
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Introduction 

The purpose of the flow sensor program is to determine if the solid state flow 
sensors can meet the target specifications for accuracy and response time and 
to characterize their overall performance (not including EMI and vibration) within 
the NASA defined operating conditions as shown in Table 1. 

The purpose of the test plan is to define the testing to be conducted to enable 
accuracy, response time data and sensor characterization to be obtained. The 
testing of the flow sensors will be divided into two categories, design of 
experiments (DOE), to obtain the sensor characterization, and one factor at a 
time testing to obtain the accuracy (linearity, hysteresis and repeatability) and 
response time data. The designed experiments will be tailored to obtain the 
maximum amount of information as inexpensively and as quickly as possible. 
The DOEs will determine sensor sensitivity factors and interaction effects for all 
of the pertinent parameters. The one factor at a time test will also tend to act as 
a verification of the DOE testing to further boost the test confidence level. 

Figure 1 shows the proposed test configuration. A data acquisition system will 
be used to both control the experimental conditions and to collect the data. 

Three DOEs will be run on the sensors. The first will be an L8 run on one > 
sensor only to carry out some preliminary screening and hence enable more 
efficient and accurate data to be subsequently collected by enabling some of the 
potential factors to be eliminated. Then a three level LI 8 DOE will be used in 
order to obtain curvature data for the sensitivity factors. Finally a second L8 will 
be run as a full factorial, again on one sensor only, to test for interactions on the 
main factors. For the LI 8 the same test conditions will be applied to both 
sensors in order to obtain replicate information which will be used to improve the 
overall accuracy of the experiment and as a confidence factor. 

The factors considered in the DOEs will be ambient pressure, ambient 
temperature, relative humidity, gas flow rate, gas temperature, gas pressure, 
gas type (three gases used in the DOEs will be oxygen, nitrogen and helium) 
and orientation about each of the X and Y axes. (Figure 2 shows the orientation 
axes definition.) These nine parameters are judged to be the only ones which 
might affect operation of the sensor. The responses to be collected and 
analyzed will be flow rate, pressure drop, power consumption, flow rate error, 
thermal conductivity error, pressure and temperature. 

The experimental factor test levels (conditions) are selected to be slightly 
different for the three experiments so that the data can be merged to give more 
complete curvature indications by having 7 points on the curve instead of only 3. 

The one factor at a time testing will consist of flow rate sweeps for each of the 
three gases in which the flow rate will be both increased from minimum to 
maximum and then decreased from maximum to minimum to obtain hysteresis 
data. Also, flow rate step response tests will be carried out for each of the three 
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gases . The one factor at a time tests will be repeated at NASA for a fourth gas 
(hydrogen) after all testing at Carleton is completed. 



Parameter 

Requirement 



Operating Fluids 

Nitrogen/Oxygen/Hydrogen/Helium 

Flow Range 

0.01 to 150 Ib/hr 

Pressure Range 


Max Pressure Drop Oxygen 



HamMa m 1 11 n im— — 

Max Pressure Drop Nitrogen 

0.5 psid @ 200 psig & 15 Ib/hr 1 



Max Pressure Drop Hydrogen 

0.5 psid @ 200 psig & 15 Ib/hr 

Max Pressure Drop Helium 

0.5 psid @ 200 psig & 15 Ib/hr 

Gas Temperature Range 

35 to 1 75 degrees F 

Ambient Pressure Range 

0 to 1 5 psia 

Ambient Temperature Range 

40 to 90 degrees F 

Accuracy 

< ± 3% of reading 

Offset at Zero Flow 

< ± 0.01 Ib/hr 

Drift 

Unit remains in calibration > 5 years 

Response Time 

< 3 seconds 

Maximum Power 

1.5 Watts® 28 VDC 

Output Voltage Range 

0 to 5 VDC 

Acceleration 

• 4.4 vertical, 4.8 lateral g’s 

EMI 

* TBD 

Volume 

• < 30 cubic inches 

Weight 

*<1.35 lbs 

Particle Impact Sensitivity 

* Oxygen safe 

Service Life 

* > 10 years 


Table 1 - Flow Sensor Requirements (* Not Required for Demonstration) 
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Measurement Accuracy for all Test Parameters 

The following list shows the available measurement accuracy for each of the 
test parameters. 

Ambient pressure 
Ambient temperature 
Relative humidity 
Gas pressure 
Gas temperature 
Flow rate 


±0.016psia 
+ TTfegree F ” 
± 0.5 % 

± 1 psig 
± 2 degrees F 
± 0.7% of reading 


2 • 3 (sin, 


Detailed Test Description For DOEs 

L8 Screening Experiment 

With nine factors to be considered, in order to obtain interaction and 
curvature data, a large L81 experiment utilizing excessive resources would be 
required. However, it is suspected that orientation about the X axis, relative 
humidity and possibly ambient pressure are less likely to affect the sensor 
performance so this L8 will be designed to find out if this is, in fact, so and hence 
enable the deletion of these factors from the main experiment. Also, flow rate, 
gas pressure and gas temperature will be included to guard against the 
possibility of interactions between these factors (which are believed to be the 
primary factors affecting sensor performance) and the factors which we wish to 
drop from the main experiment. 

This L8 will be run with only nitrogen and only one sensor because it is 
believed that sensor performance with one gas and one sensor will be 
characteristic and show the same trends with any gas or any sensor. 

Figure 3 shows the actual experimental treatment combinations and levels. 

LI 8 Main Experiment 

Assuming that the L8 screening experiment verifies the expectation that 
orientation about the X axis, relative humidity and ambient pressure can be 
eliminated then the LI 8 main experiment will obtain curvature and sensitivity 
data for gas type, flow rate, gas pressure, gas temperature, ambient temperature 
and orientation about the Y axis for each of the two sensors. 

Figure 4 shows the actual experimental treatment combinations and levels. 
This experiment can accept up to 7 three level parameters and one 2 level 
parameter so that with only 6 parameters actually utilized a measure of 
experimental error can be obtained because the experiment is not saturated. 
However, because any interactions between any pair of parameters in columns 3 
through 6 are partially confounded with columns 3 through 8 there is a risk that if 


file : FLOW .DAT 

Title: Solid State Flow Sensor Screening Test 
Comment: Preliminary Sensitivity Identification 
Re^nonse: Flow & Err, Pwr, Delta P & Thermal Err 



n Listing: 


TC# 

12 3 4 



Gas 

Flow 

Press 

Gas Temp 

1 

Oxygen 

1 

40 

30 

2 

Oxygen 

10 

475 

105 

3 

Oxygen 

100 

910 

180 

4 

Nitrogen 

1 

40 

105 

5 

Nitrogen 

10 

475 

180 

6 

Nitrogen 

100 

910 

30 

7 

Helium 

1 

475 

30 

8 

Helium 

10 

910 

105 

9 

Helium 

100 

40 

180 

10 

Oxygen 

1 

910 

180 

11 

Oxygen 

10 

40 

30 

12 

Oxygen 

100 

475 

105 


^ Nitrogen 

1 

475 

180 

fl 

m Nitrogen 

10 

910 

30 

n 

W Nitrogen 

100 

40 

105 

16 

Helium 

1 

910 

105 

17 

Helium 

10 

40 

180 

18 

Helium 

100 

475 

30 


Gas 

Flow 

Press 

Gas Temp 


(L*- 



8 


Y Attitude 


Up 

Up 

Up 

Up 

Up 

Up 

Up 

Up 

Up 

Down 

Down 

Down 

Down 

Down 

Down 

Down 

Down 

Down 


Y Attitude 


Figure 4 



the interactions exist and are large then experimental error will be confused with 
them. Thus, a full factorial, follow on experiment will be run to identify the 
interactions. 

L8 Full Factorial Experiment 

The assumption is made that if interactions exist they are most likely between 
flow rate and gas temperature and flow rate and gas pressure. Thus with three 
factors an L8 is required for no fractionalization to occur. If the LI 8 shows any 
possible interactions between orientation and any parameter then this 
experiment may have to be expanded to an LI 6. Figure 5 shows the actual 
experimental treatment combinations and levels. Like the L8 screening 
experiment this experiment only needs to be run for nitrogen and one sensor 
because it again is assumed that trends for one gas and one sensor are 
characteristic for all gases. 

One Factor At A Time Tests 

These tests will utilize the same test setup as shown in figure 1 and will be 
run with ambient temperature at normal factory room temperature of 
approximately 72 degrees F ± 7 degrees and with orientation normal with the X 
axis horizontal (parallel to the factory floor). Each flow sensor will be subjected 
to the same tests. The tests will be repeated for each of the three gases 
nitrogen, oxygen and helium. Each test will consist of a flow sweep from zero 
flow up to maximum flow for the particular gas and back to zero while attempting 
to maintain gas temperature and gas pressure relatively constant. The test will 
be run at each of four conditions of gas temperature and pressure comprising all 
combinations of high and low. Thus high gas pressure and temperature, low gas 
pressure and temperature, high gas pressure and low gas temperature and 
finally, low gas pressure and high gas temperature. This will then require a total 
of 24 flow sweeps to complete all tests. This test data will be analyzed for 
linearity, error and hysteresis. 

The test setup will then be modified to place a quick opening ball valve 
between the flow sensor and the pressure regulator. The data acquisition data 
collection rate will be set to maximum and the number of parameters collected 
will be reduced to flow rate and gas pressure only to facilitate data collection. 
The ball valve will be rapidly opened and then closed to simulate rising and 
falling step inputs to the sensor. The collected data will be analyzed to estimate 
the sensor response time. This series of tests will be repeated at an average gas 
temperature only but otherwise will cover the same combinations and conditions 
as the flow sweep tests and so will require 12 tests to complete the sequence. 
The volume between the ball valve and the sensor and between the sensor and 
the handvalve will be measured and hence by knowing the pressure rise and fall 
rates of these volumes the theoretical flow rates versus time can be computed 
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and compared with the actual and by this means obtaining a measure of sensor 
response time. 




FLOW SENSOR TEST PROCEDURE 




Required Test Condition Tolerances for all Experiments and Tests 


Ambient pressure 
Ambient-temperature- 
Relative humidity 
Gas pressure 
Gas temperature 
Flow rate 
Gas purity 


±0.25 psia 
±2 dogrooo - F 
±5% 

± 8 psig 
± 5 degrees F 
± 1 % of setting 
± 0.1 % 


RA . 
?/?/u 


Measurement Accuracy for all Test Parameters 


Ambient pressure 
Ambient temperature 
Relative humidity 
Gas pressure 
Gas temperature 
Flow rate 

L8 Screening Experiment 


±0. 016 psia 

±^3egreiT~ ^ , 

±0.5% 

± 1 psig 
± 2 degrees F 
± 0.7% of reading 


1 . Use only the Honeywell unit A for this experiment. 

2. Setup the test equipment in the configuration shown in figure 1 . 

3. To ensure thermal stability, allow the Honeywell unit and the reference 
flow sensor to warm up for 10 minutes after power up before any readings are 
taken. 

4. Use a multi-turn handvalve with a maximum orifice capable of flowing 
1000 LPM at 40 psig. 

5. Set the DAC to collect data once every 10 seconds. Ensure that the DAC 
data is saved to disk, ensure that all actual test condition data is recorded as 
well as the electrical current to the unit under test and the following outputs from 
the unit; flow rate, pressure, temperature and thermal conductivity. 

6. Use nitrogen as the supply gas. 

7. For the reference flow sensor use the laminar flow elements for 5 and 50 
Ib/hr ranges (refer to table 1 for the correct elements). For the first test in the 
experiment use the 5 Ib/hr element. 

8. Refer to figure 3 for the test conditions required for this screening 
experiment. Note that the screening tests may be conducted in whichever order 
is most economical. 
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HONEYWELL FLOW SENSOR TEST SETUP 










METER NUMBER 






RANGE (cfm) 

RANGE (Ib/hr) 

RANGE (Ib/hr) 

RANGE (Ib/hr) I 


LOW 

HIGH 

LOW 

HIGH 

LOW 

HIGH 

LOW 

HIGH | 





bhhhhh 

1 

0.002 


0.0088 

0.0352 

0.0100 

0.0403 

0.0012 


2 

0.008 


0.0352 

0.1418 

0.0403 

0.1621 

0.0050 

0.0202 

3 

0.033 


0.1418 

0.5710 

0.1621 

0.652 

0.0202 

0.081 

4 

0.131 

0.5287 

0.5710 

2.2990 

0.6525 

2.627 

0.0815 

0.328 

5 

0.529 

2.1286 

2.2988 

9.2550 

2.6268 

10.5754 

0.3281 

1.321 

6 

mam 

8.5695 

9.2545 

37.2580 

10.5754 

42.575 

T'ltri 

5.317 

7 

8.57 

34.5000 

37.2579 

150.0000 

42.5754 

171.4000 

5.3170 

21.410 


Table 1 - Flow Meter Ranges 



9. Establish on the Honeywell unit A the location of the sensing element (the 
T018 can on the side of the venturi) and arrange the unit in the PTH chamber so 
that the T018 can is pointing in the up direction per figure 6 for all tests where 
the X attitude is required to be +90 degrees and in the down direction for tests 
where the attitude is -90 degrees. 

1 0. Adjust the pressure regulator for 1 00 psig and adjust the hand valve for 5 
Ib/hr on the reference flow sensor. Re-adjust the pressure to compensate for any 
droop. 

1 1 . Adjust the PTH chamber to the correct conditions, ambient pressure 1 4.5 
psia and relative humidity 85 %. Maintain the chamber temperature at 70 
degrees F at all times. 

1 2. Adjust the temperature of the gas heater/chiller until the gas temperature 
is 50 degrees F and recheck the flow and pressure. Allow the unit to stabilize for 
2 minutes after all test conditions are stabilized. 
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1 3. Adjust the test conditions for the next test and repeat the above 

procedure, steps 9 through 12 being careful to always use the correct laminar 
flow element. 


LI 8 Main Experiment 

1 . Use both the Honeywell units A and B in turn for this experiment. 

2. Setup the test equipment in the configuration shown in figure 1 . 

3. To ensure thermal stability, allow the Honeywell unit and the reference 
flow sensor to warm up for 1 0 minutes after power up before any readings are 
taken. 

4. Use a multi-turn handvalve with a maximum orifice capable of flowing 
1000 LPM. 

5. Set the DAC to collect data once every 1 0 seconds. Ensure that the DAC 
data is saved to disk, ensure that all actual test condition data is recorded as 
well as the electrical current to the unit under test and the following outputs from 
the unit; flow rate, pressure, temperature and thermal conductivity. 

6. For the reference flow sensor use the laminar flow elements for 1 ,1 0 and 
100 Ib/hr ranges (refer to table 1 for the correct elements). For the first test in the 
experiment use the 1 Ib/hr element. 

7. Refer to figure 4 for the test conditions required for this main experiment. 

Note that the individual tests may be conducted in whichever order is most 
economical. 

8. Arrange the plumbing in the chamber so that the Honeywell units can be 
installed with the flow direction pointing up or down as required for the test 
configuration. See figure 7 for orientation. Initially install in the up direction for 
the first test. 

9. Connect the oxygen (or nitrogen or helium) gas supply to the pressure 
regulator. 

1 0. Adjust the pressure regulator for 40 psig. 

1 1 . Adjust the hand valve for 1 Ib/hr on the reference flow sensor. Re-adjust 
the pressure to compensate for any droop. 

1 2. Adjust the PTH chamber to the correct condition, temperature aL5S' 3c \ . 

degrees F. Maintain the ambient pressure at approximately 1 4.7 psia and ^ / 

relative humidity at approximately 30 % at all times, try and keep constant. ' 

1 3. Adjust the temperature of the gas heater/chiller until the gas temperature 
is 30 degrees F and recheck the flow and pressure. Allow the unit to stabilize for 
2 minutes after all test conditions are stabilized. 

1 4. Adjust the test conditions for the next test and repeat the above 
procedure, steps 8 through 13 being careful to always use the correct laminar 
flow element. 



1 . Use only one of the Honeywell units A or B for this experiment. (Which 
one will be decided after main experiment test results are obtained.) 

2. Setup the test equipment in the configuration shown in figure 1 . 

3. To ensure thermal stability, allow the Honeywell unit and the reference 
flow sensor to warm up for 10 minutes after power up before any readings are 
taken. 

4. Use a multi-turn handvalve with a maximum orifice capable of flowing 
1000 LPM at 40 psig. 

5. Set the DAC to collect data once every 1 0 seconds. Ensure that the DAC 
data is saved to disk, ensure that all actual test condition data is recorded as 
well as the electrical current to the unit under test and the following outputs from 
the unit; flow rate, pressure, temperature and thermal conductivity. 

6. Use nitrogen as the supply gas. 

7. For the reference flow sensor use the laminar flow elements for 2 and 80 
Ib/hr ranges (refer to table 1 for the correct elements). For the first test in the 
experiment use the 2 Ib/hr element. 

8. Refer to figure 5 for the test conditions required for the full factorial 
experiment. Note that the tests may be conducted in whichever order is most 
economical. 

9. Plumb the unit under test so that the flow direction is horizontal and the 
sensing element is also horizontal (Y axis orientation). 

1 0. Adjust the pressure regulator for 70 psig and adjust the hand valve for 2 
Ib/hr on the reference flow sensor. Re-adjust the pressure to compensate for any 
droop. 

1 1 . Adjust the PTH chamber to the correct conditions, ambient pressure 1 4.7 
psia, temperature at 70 degrees F and relative humidity 30 %. Maintain the 
chamber at these conditions at all times. 

1 2. Adjust the temperature of the gas heater/chiller until the gas temperature 
is 40 degrees F and recheck the flow and pressure. Allow the unit to stabilize for 
2 minutes after all test conditions are stabilized. 

1 3. Adjust the test conditions for the next test and repeat the above 
procedure, steps 10 through 12 being careful to always use the correct laminar 
flow element. 
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File: FLOW 2 .DAT 

Title: Solid State Flow Sensor Screening Test 

Comment: Identification of Interactions 

Response: Flow & Err, Pwr, Delta P & Thermal Err 

Listing: 



TC# 


1 


2 


3 


Flow Rate 

1 : 2 
2: 80 
3 : 2 

4 : 80 

5: 2 

6: 80 
7: 2 

8 : 80 

Flow Rate 


Pressure 

Gas 

70 

40 

70 

40 

850 

40 

850 

40 

70 

165 

70 

165 

850 

165 

850 

165 

Pressure 

Gas 


Temp 


Temp 


Figure 5 


UP 



Figure 6 X Axis Orientation 


UP 

FLOW 

DIRECTION 



T 

DOWN 

FLOW 

DIRECTION 



Figure 7 Y Axis Orientation 


One Factor At A Time Tests 


Linearity, Accuracy and Hysteresis 

1 . Use both the Honeywell units A and B in turn for these tests. 

2. Setup the test equipment in the configuration shown in figure 1 . 

3. To ensure thermal stability, allow the Honeywell unit and the reference 
flow sensor to warm up for 1 0 minutes after power up before any readings are 
taken. 

4. Use a multi-turn handvalve with a maximum orifice capable of flowing 
1000 LPM at 40 psig. 

5. Set the DAC to collect data once every 1 second. Ensure that the DAC 
data is saved to disk, ensure that all actual test condition data is recorded as 
well as the electrical current to the unit under test and the following outputs from 
the unit; flow rate, pressure, temperature and thermal conductivity. 





6. Plumb the unit under test so that the flow direction is horizontal and the 
sensing element is also horizontal (Y axis orientation- see figure 2). 

7. For the reference flow sensor use the appropriate laminar flow elements 
for the flow range (refer to table 1 for the correct elements). Note that this test 
will probably require the use of all nine elements in turn. 

8. Connect the nitrogen, oxygen and helium gas supply in turn as required to 
the pressure regulator. 

9. Adjust the pressure regulator for 40 psig with the hand valve closed. 

1 0. Adjust the PTH chamber to the correct conditions, ambient pressure 1 4.7 
psia, temperature at 72 degrees F and relative humidity 30 %. Maintain the 
chamber at these conditions at all times. 

1 1 . Open the handvalve to obtain some flow and adjust the temperature of the 
gas heater/chiller until the gas temperature is 30 degrees F. Allow the unit to 
stabilize for 2 minutes after all test conditions are stabilized. 

1 2. Close the handvalve and start the data recording. Slowly open the 
handvalve and increase the flow. Keep checking the pressure is constant. When 
the first laminar flow element has reached its maximum stop the handvalve 
motion and change elements. Continue increasing flow. When full flow of 150 
Ib/hr is reached reverse the process and start slowly reducing flow until zero is 
reached. 

13. Test unit B in the same way. 

1 4. Set the next test conditions for pressure and gas temperature (refer to 
table 2) and repeat the flow sweep. 

1 5. After all table 2 tests are complete for nitrogen for both units A and B 
connect oxygen as the supply gas and repeat the table 2 tests. 

1 6. Finally, repeat again for helium as the supply gas. 


Pressure 

Temperature 

40 

30 

40 

180 

900 

30 

900 

180 


Table 2 


Response Time 

1 . Modify the test setup to introduce a quick opening ball valve between the 
pressure regulator and the Honeywell unit. 



2. Modify the DAC setup to record only the Honeywell unit parameters and 
the reference flow sensor parameters. Increase the DAC sample rate to its 
maximum, or 20 times per second, which ever is smaller. 

3. Set the gas temperature to approximately 1 05 degrees F. 

4. At a gas pressure of 40 psig for nitrogen and with Honeywell unit A 
quickly open, hold for 5 seconds and then quickly close the ball valve. 

5. Repeat step 4 at a pressure of 900 psig. 

6. Repeat both steps 4 and 5 (pressures of 40 and 900 psig respectively) 
for unit B. 

7. Change the supply gas to oxygen and repeat the four tests. 

8. Change the supply gas to helium and again repeat the four tests. 
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Abstract 


The Space Shuttle Orbiter Atmospheric 

« n o\/jtalization Pressure Control System (ARPCS) and 
Fuel Cell System (FCS) use a hot wire anemometer 
of gas mass flow sensor for flow measurement. In 
ARPCS oxygen and nitrogen mass flows are 
measured and in the FCS oxygen and hydrogen mass 
flows are measured. The existing flow sensors suffer 
from certain accuracy limitations and potential failure 
modes. A new type of commercially developed solid 
state micro-machined silicon gas mass flow sensor 
developed by Honeywell was adapted to allow the 
technology to be assessed for the application. 

A demonstration test program has been 
conducted to evaluate the performance characteristics 
of the new sensor for space system applications and 
environments. The testing was sponsored by the 
National Aeronautics and Space Administration (NASA) 
at the Johnson Space Center (JSC). The primary testing 
was conducted at Carleton Technologies Inc. using 
nitrogen, oxygen and helium (as a surrogate for 
hydrogen). Hydrogen testing will be conducted at JSC. 
The objective of the tests was to determine if the new 
flow sensing technology would benefit the Shuttle and 
Space Station programs. 

The new “smart” sensor has many advantages 
over the existing technology. These include low power 
consumption, very low drift rate, measurement 
capability over a 15000:1 flow range (compared to less 
ip 1000:1 for the existing technology), high accuracy 
the ability to be used with different gas types 



without modification or recalibration. The sensor has 
built in temperature and gas composition compensation, 
based upon heat transfer measurements. 

The paper discusses the sensor operating 
principle, the test methodology, the sensor operational 
requirements and the test results. 

Introduction 


The Space Shuttle Orbiter Atmospheric 
Revitalization Pressure Control System (ARPCS) and 
the Fuel Cell System (FCS) presently use a hot wire 
anemometer type of gas mass flow sensor for flow 
measurement. In the ARPCS oxygen and nitrogen mass 
flow are measured and in the FCS oxygen and hydrogen 
mass flow are measured. The measurement method 
utilized consists of sensing the amount of heat 
transferred from a heater to a sensing element by the 
flow of gas. The existing sensors suffer from three main 
deficiencies which affect their accuracy and have 
necessitated a search for a superior product. The 
problems experienced by the existing sensors are 
orientation sensitivity, limited range and long term drift. 
The orientation effect is caused by thermal micro- 
currents which are sensitive to the gravitational vector 
and hence to position. In addition, the existing unit offers 
an accuracy of 3% of full scale over a range of 0 to 5 
Ibs/hour while the goal is to achieve 3% of reading (with 
a minimum resolution of ± 1% of reading) over a range 
of 0.01 to 150 Ibs/hour. 

The flow sensing applications on the Space 
Shuttle present very severe requirements for almost any 
presently available technology. A recently developed 
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Figure t - Flow sensor type 1 


micro-machined solid state silicon technology appears 
*o have the potential to meet these requirements. Due to 

« |mall size, low heat dissipation and its inert material 
^position (Pt and Si 3 N 4 ) it was expected that it would 
nificantly reduce, or eliminate, both the orientation 
and long term drift errors. To verify this a program to 
demonstrate the ability of the new design to meet the 
Shuttle ARPCS and FCS accuracy, pressure drop, 
power consumption and response time requirements 
was conducted. In this program two types of existing 
“smart” sensor structures 0 ' 7 ' were adapted and tested 
over a wide environmental range. Originally testing of 


only one type was envisaged but disappointing 
performance in certain areas led to the inclusion of the 
second type. 

Type 1 consisted of one flow microsensor 
mounted at the throat of a Venturi nozzle 0 ' 3 ’ while type 2 
had the same sensor mounted in a bypass* 4 ' 5 ’. The latter 
approach allows the use of a commercial sensor 
package and ensures that the microsensor: (1) is only 
exposed to flows driven by a multi-tap, averaged 
(piezometric sampling) pressure drop across a laminar 
flow restriction (honeycomb): and (2) is factory- 
assembled reproducibly in a well-defined flow 
microchannel. This approach also permits a wide 
enough flow range that the special high-flow 
measurement method (based on the measurement of 
heater power dissipation) could be eliminated, reducing 
^«|jdware, software and calibration costs. Figures 1 and 
^^Bpow cross-sectional sketches of the respective flow 
^^reor configurations. 


Following the successful completion of the proof 
of concept program (with satisfactory accuracy 
performance measurements) the next stage would be a 
full development and qualification program during which 
other issues would be resolved, such as EMI and 
vibration. 

This paper describes briefly the principle of 
sensor operation (a full description is available in past 
papers’ 1 ' 7 ’ ), the requirements imposed on a flow sensor 
in the Space Shuttle environment, the strategy adopted 
for a short but thorough test of the concept, the test 
results obtained for both sensor configurations and 



Figure 2 - Flow sensor type 2 

finally, conclusions and recommendations for further 
work. 

Flow Sensor Principle of Operation 

The silicon based microbridge sensor uses a hot 
film element and a differential sensing approach that is 
more sensitive than existing hot wire anemometer 
designs. The “smart” microstructure sensor is micro- 
processor controlled and consists of two independent 
thermal sensing units. One sensing unit is placed 
upstream of a Venturi metering channel (figure 3), and 
the other is placed either at the throat of the Venturi for 
Type 1 (figurel) or in a bypass for Type 2 (figure 2) . 
The upstream sensing unit, or the property sensor, is 
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Gas Property Sensor 



Figure 3 - Venturi flow channel 

used to evaluate the gas properties (i.e. the 
temperature, the thermal conductivity and the specific 
heat of the gas). This allows the sensor output to be 
automatically compensated for different gas (1) types. 
The property sensor is recessed to avoid direct flow 
impingement, which would disturb its measurement 
function. The structure of the property sensor is identical 
to that of the flow sensor. 

The flow sensing unit is used to measure the 
uncompensated flow. Each flow sensor consists of a 
heater located in the middle of two sensing elements 
which are positioned upstream and downstream with 
respect to the direction of gas flow (figure 4). The sensor 
can be operated in two modes, one of which gives best 
accuracy at I6w flows and one which gives best 
accuracy at high flows. 

The first mode is that in which the differential signal 
between the upstream and downstream sensing 
^ments is used. In the second mode the total heat 
[^ipation of the heating element is used. The 
l^oprocessor merges the two flow outputs to generate 
a smooth changeover from one to the other as a 
function of flow. 



Temperature Sensor 



Figure 4 - Chip layout 


The micro-processor also receives outputs from 
the property and temperature sensors and determines 
the values of k (thermal conductivity), c p (specific heat) 
and temperature. It also receives the uncompensated 
signals from the flow sensor. The mass flow rate of the 
gas is then computed by correcting the flow signal for 
pressure, temperature, c p and k. The sensor was 
configured with a 0 to 5 volt output to match the signals 
levels from the existing Shuttle sensor. However, 
because the existing sensor maximum range is only 5 
Ib/hr and the microbridge sensor maximum range is 150 
Ib/hr it was advantageous to convert the final micro- 
processor output to a voltage logarithmically 
proportional to flow so as to maintain the resolution and 
hence the required accuracy for low flows. 

Sensor Characteristic Objectives 

Listed below in Table 1 are the performance 
goals for the sensor. Table 2 specifies the 

environmental operating regime and in Table 3 the 
physical characteristic requirements are outlined. 

Table 1 - Flow Sensor Performance Goals 


Parameter 

Requirement 



Flow Ranqe 

0.01 to 150 Ib/hr 

Accuracy 

< ± 3% of readinq 

Offset at Zero Flow 

<±0.01 Ib/hr ■ 

Drift 

Unit remains in calibration > 5 
years 

Max Pressure Drop Oxyqen 

0.5 psid @ 900 psiq & 15 Ib/hr 


2.0 psid @ 900 psiq & 75 Ib/hr 

Max Pressure Drop Nitroqen 

0.5 psid @ 200 psiq & 15 Ib/hr 


28.0 psid @ 200 psiq & 150 Ib/hr 

Max Pressure Drop Hydroqen 

0.5 psid @ 200 psiq & 15 ib/hr 

Max Pressure Drop Helium 

0.5 psid @ 200 psiq & 15 Ib/hr 

Response Time 

< 3 seconds (T90) 

Maximum Power 

1.5 Watts @ 28 VDC 

EMI 

TBD 

Service Life 

> 1 0 years 


Table 2 - Flow Sensor Environmental Operating 
Regime 


Parameter 

Requirement 



Operatinq Fluids 

Nitroqen Oxyqen Hydroqen Helium 

Gas Temperature Ranqe 

35 to 175 deqrees F 

Ambient Pressure Ranqe 

0 to 1 5 psia 

Ambient Temperature Ranqe 

40 to 90 deqrees F 

Acceleration 

4.4 vertical. 4.8 lateral q's 


It was decided that weight, volume, particle 
impact sensitivity, acceleration and EMI were outside 
the scope of proof of concept testing and were not 
included in the test program. 
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Table 3 - Flow Sensor Physical Characteristic 
Requirements 



The purpose of the flow sensor program was to 
determine if the solid state flow sensors were able to 
meet the target specifications for accuracy, pressure 
drop, power consumption and response time and to 
characterize the overall performance (not including EMI 
and vibration) within the NASA defined environmental 
regime shown in Table 2. 



OAC . DATA ACQUSTION 
VaHANOVALVt 

Figure 5 - Test Setup 


The testing of the flow sensors was divided into 
two categories, Taguchi style statistical experimental 
designs (SED) to characterize the sensors, and “one 
factor at a time” testing to obtain the accuracy (linearity, 
hysteresis and repeatability) and response time data. 
The designed experiments were tailored to obtain the 
maximum amount of information as inexpensively and 
as quickly as possible. The SEDs were used to 
determine sensor sensitivity factors and interaction 
effects for all of the pertinent parameters. The one 
factor at a time tests also validated the SED testing to 
further boost the test confidence level. 

Statistical experimental designs are used to 
determine variation due to predetermined parameters 
under study. These parameters are controlled in a 
manner such that a statistical evaluation will pinpoint 
which parameters or combinations of parameters are 
causing certain effects. An experiment matrix should be 
designed so there are enough columns (one for each) to 
accomodate the parameters under study, interactions 

P " 1 ^ results of combinations of parameters) and 
erimental error. If the matrix is too small (saturated) 
interactions and experimental error will be confused 
with individual parameters. 


Three SEDs were planned. The first was a 
preliminary screening test, using one sensor, to identify 
unimportant factors which could then be eliminated from 
subsequent testing making this testing more efficient. A 
three level SED was used to obtain curvature data for 
the sensitivity factors. Finally one sensor was subjected 
to a full factorial SED to test for interactions of the main 
factors. The curvature experiment was performed twice, 
once for each sensor configuration. The test conditions 
were kept constant to obtain replicate information which 
was used to improve the overall accuracy of the 
experiment and as a confidence factor. 

The factors considered in the SEDs were 
ambient pressure, ambient temperature, relative 
humidity, gas flow rate, gas temperature, gas pressure, 
gas type (three gases used in the SEDs were oxygen, 
nitrogen and helium) and orientation about the X and Y 
axes. These nine parameters are judged to be the only 
ones which might affect the accuracy of the sensor in 
the test scenario defined. The responses to be collected 
and analyzed were flow rate, pressure drop, power 
consumption and flow rate error. 

The ranges of the experimental factor test levels 
(conditions) were selected for the three experiments to 
be similar, but not identical, so that the data could be 
merged to give more complete curvature indications. 

The “one factor at a time” testing consisted of 
flow rate sweeps at all combinations of pressure- and 
temperature limits for each of the three gases. For these 
tests the flow rate was increased from minimum to 
maximum and then decreased from maximum to 
minimum to obtain hysteresis data. Also, flow rate step 
response tests were carried out for each of the three 
gases . The one factor at a time tests will be repeated at 
Johnson Space Center for a fourth gas (hydrogen) after 
all testing at Carleton is completed. 

Figure 5 shows the test configuration. A data 
acquisition system was used to both control the 
experimental conditions and to collect the data. 

Because the result of the test program was 
intended to assess the accuracy of the flow sensor, the 
accuracy of the test equipment used was critical. This 
was especially true for the flow standard used. A series 
of seven laminar flow elements were used to cover the 
flow range from 0.01 to 150 pounds per hour for oxygen, 
nitrogen and helium. Laminar flow elements were 
chosen because of their ability to respond as fast as the 
smart sensor. The particular elements were chosen to 
give an accuracy of 3 to 4 times better than the 
specification flow accuracy requirement. Table 4 shows 


Table 4 - Measurement Accuracy 


Ambient pressure 

± 0.016 psia 

Ambient temperature 

± 1 deqree F 

Relative humidity 

±0.5% 

Gas pressure 

± 1 psiq 

Gas temperature 

: ± 2 deqrees F 

Flow rate 

± 0.7% of reading 


the measurement accuracy for each of the test 
parameters. 
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The LI 8 three level main experiment was 
designed to obtain curvature and sensitivity data for gas 
type, flow rate, gas pressure, gas temperature, ambient 
temperature and orientation about the Y axis for each of 
a two sensors. 



Table 5 - LI 8 Experimental Conditions 


Gas 

Flow 

Pressure 

Temp 

Attitude 

Oxygen 

1 

40 

30 



Oxygen 

10 

475 

105 

Up 

Oxygen 

100 

910 

180 

Up 


1 

40 

105 

_ue 

Nitrogen 

10 

475 

180 

Up 

Nitrogen 

100 

910 

30 

_ye 

Helium 

0.1 

475 

30 

.ye 

Hefium 

1 

910 

105 

Up 

Helium 

10 

40 

180 

Up 

Oxygen 

1 

910 

180 

Down 

Oxygen 

10 

! 40 

30 

Down 

Oxygen 

100 

i 475 

105 

Down 

Nitrogen 

1 

i 475 

180 

Down 

Nitrogen 

10 

910 

30 

Down 

Nitrogen 

100 

40 

105 

Down 

Helium 

0.1 

910 

105 

Down 

Helium 

1 

40 

180 

Down 

Helium 

10 

475 

30 

Down 


Table 5 shows the actual experimental 
treatment combinations and levels. This experiment can 
accept up to 7 three level parameters and one 2 level 
parameter so that with only 5 parameters actually 
utilized, a measure of experimental error can be 
htained because the experiment is not saturated. 

f fct/ever, because any interaction between pressure 
■ temperature is partially confused with Y attitude and 
e spare columns used to collect experimental error 
there is a risk that if this interaction exists then 
experimental error will be confused with it. Thus, full 
factorial, follow on experiments were planned to identify 
interactions. 


NITROGEN @ 200 PSIG AND 70 DEGREES F 


1000 



Figure 6 - Unit 1 Preliminary Flow Test 

on the sensor output as expected. The test results are 
shown in Table 6. 

The ANalysis Of VAriance (ANOVA) of the flow 
signal showed that 98.43% of the variation in flow was 
due to the flow changes within the experiment which 
indicates, not unexpectedly, that as a flow sensor the 
unit is very good. It also shows that the other factors 
within the experiment, humidity, gas pressure, gas 
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Figure 7 - Unit 1 Flow Error 


Test Results Unit 1 

Prior to commencing the main test program a 
preliminary flow test, using nitrogen at room 
temperature, was conducted to verify correct operation 
of the sensor and the test equipment. The results are 
shown in Figures 6 and 7. These show that the sensor 
gives an output with a linear trend over the required flow 
range but with a maximum deviation of 21% of reading. 

It was determined that, apart from a narrow flow 
range, the deviation was significantly constant and thus 
could be attributed to a correctable calibration error. 
Therefore the test program was continued knowing that 
absolute errors would be unacceptable but in the 
knowledge that average error could be improved by 
more accurate calibration. 

The first structured experiment (L8), which was 
designed to verify that humidity, ambient air pressure 
X-axis orientation did not influence performance, 
successful in that analysis of the test results 
onstrated that these factors had a negligible effect 


temperature, ambient pressure and X-axis orientation 
produced output variation of only 1.42% and so had 
almost no effect on the flow signal. The experimental 
error was only 0.15% of the variation within the 
experiment 

Table 6 - Screening Test Results 


RUN 

NUMBER 

ACTUAL 

FLOW 

RS232 

FLOW 

ERROR % OF 
READING 

1 

5.06 

5.307 

4.88 

2 

50.6 

54.2 

7.11 

3 

5.02 

5.375 

7.07 

4 

50.65 

54.2 

701 

5 

5.015 

5.279 

5.26 

6 

50.14 

57.53 

14.74 

7 

5.074 

5.245 

3.37 

8 

49.36 

65.7 

33.1 


showing that the conditions were accurately set and that 
there were no significant effects from outside sources. 
An ANOVA was also performed on the flow error (as a 
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percentage of reading) and showed that 31.8% of the 
variation in error was due to flow, 6.4% due to gas 
pressure, 17.1% due to gas temperature, 6.0% due to 
ambient air pressure, 25.5% due to humidity and 3.8% 
•'••e to X-axis orientation with 9.4% experimental error. 

* was hypothesized that the apparent effects, 
Icularly the large and unexpected effect of humidity 
the error were due to the use of a saturated 
experiment and to the difficulty in maintaining consistent 
gas and unit temperatures. To verify this a full factorial 
L4 experiment was run with flow and humidity as the 
only factors. The results are shown in Table 7. The 
ANOVA showed that the flow changes accounted for 
99.9% of the variation in flow and humidity for only 
0.05%. The effect of flow error due to humidity in this 
experiment was not statistically significant. 


Table 7 - Humidity Test 


FLOW 

LB/HR 

HUMIDITY 
% RH 

ACTUAL 

FLOW 

RS232 

FLOW 

ERROR % OF 
READING 

5 

25 

5.02 

5.355 

6.67 

50 

25 

50.79 

54.07 

6.46 

5 

85 

5.056 

5.438 

7.56 

50 

85 

50.21 

52.05 

3.66 


The results for the LI 8 experiment are shown in 
Table 8. The flow signal ANOVA showed that 99.43% of 
the flow variation was due to flow changes within the 
experiment with no statistically significant effects from 
other sources. The flow error ANOVA showed no 
statistically significant effects, and a residual of 55.6% 
ilich indicates that there was more variation in flow 
r due to experimental error and outside factors than 
e was due to the factors within the experiment. Of 
the five factors within the experiment, the gas pressure 
and the flow level, were the main contributors to flow 




Table 8 - LI 8 Test Results 


RUN 

NUMBER 

ACTUAL 

FLOW 

RS232 

FLOW 

ERROR % OF 
READING 

t 

t .048 

1.1326 

8.07 

2 

9.9807 

10.825 

8.46 

3 

101.75 

167.5 

64.62 

4 

1.0 

1.078 

7.8 

5 

10.01 

13.725 

37.11 

6 

100.55 

105.03 

4.46 

7 

0.1067 

0.1235 

15.75 

8 

1.0023 

1.043 

4.06 

9 

9.998 

9.115 

-8.83 

10 

1.011 

1.025 

1.38 

11 

10.066 

11.118 

10.45 

12 

101.18 

140.48 

38.84 

13 

1.002 

1.0972 

9.5 

14 

9.997 

11.235 

12.38 

15 

100.64 

108.32 

7.63 

16 

0.10768 

0.1222 

13.48 

17 

0.978 

1.0261 

4.92 

18 

10.119 

15.658 

54.74 


error. The main effects plots (Figure 8 on the next 
e) show that the sensor output is significantly linear 
flow and virtually not affected by the other factors. 



The main experiment was found to have been 
poorly designed because the flow levels selected were 
not equally spaced and because the interactions were 
uniformly distributed across all columns which gave the 
indication of many interactions. The lack of equal flow 
level spacing was overcome by analysis of the 
logarithmic flow value. A series of five “two-level” 
experiments was subsequently necessary in order to 
establish which interactions were real and which were 
spurious. The interactions indicated by the main 
experiment were between gas pressure and gas 
temperature, between gas temperature and Y attitude, 
between gas pressure and Y attitude, between gas 
pressure and gas type and between gas temperature 
and gas type. Subsequent full factorial experimentation 
using L4 SED’s showed that the only real interactions 
are between gas pressure and gas temperature and 
between gas pressure and Y attitude. The results for 
these two L4 experiments are shown in Tables 9 and 10 
and the interaction plots in Figures 10 and 12. Both of 
the interactions were relatively weak being about 1.4% 
and 4.5% of flow reading respectively but together made 
it almost impossible to achieve the target 3% of reading 
requirement. 

The gas pressure and Y attitude interaction was 
caused by the gravitational effect previously discussed 
in the introduction and was identified by prior testing. 
The L4 results showed that at low pressure (1 00 psi) the 
gravitational effect was almost zero but at high pressure 
(900 psi) was quite large such that when gas was flowing 
up through the unit the measured flow rate was 
approximately 4.5% higher than when gas was flowing 
down through the unit. The gas pressure and gas 
temperature interaction was previously unknown and 
showed what appeared to be a gas density related effect 
such that at low temperature and high pressure (when 
gas density was a maximum) the unit measured flow 
rate output was about 1.4% higher than when the 
pressure was low or the temperature was high. 

Table 9 - Pressure/Temperature Interaction 


GAS 

TEMP 

GAS 

PRESSURE 

ACTUAL 

FLOW 

CORR 

RS232 

ERROR % OF 
READING 

40 

70 

1.99 

2.193 

9.65 

165 

70 

1.989 

2.184 

9.2 

40 

850 

2.024 

2.2243 

11.22 

165 

850 

1 .9928 

2.1798 

8.99 


Table 10 - Y Attitude/Pressure Interaction 


GAS 

PRESSURE 

Y 

ATTITUDE 

ACTUAL 

FLOW 

CORR 

RS232 

ERROR % OF 
READING 

100 

UP 

2.0166 

2.2037 

10.19 

850 

UP 

2.014 

2.2502 

12.51 

100 

DOWN 

2.008 

2.1982 

9.91 

850 

DOWN 

2.006 

2.1456 

7.28 


6 






Figure 8 - Unit 1 Main Effects Plot 


Plots of the smart sensor flow signal output 
versus the actual measured flow for each of the seven 
laminar flow elements (LFE’s) necessary to cover the 
flow range were obtained. Figure 9 is a concatenation of 
the seven individual tests and shows the linearity and 
hysteresis performance for one particular gas type, 
pressure and temperature combination. This plot is 
representative of the other data collected. 

Figure 11 shows the error as a percentage of 
reading. These graphs show clearly the problems 
associated with measuring flow rate over such a wide 
range and problems with this particular unit. 





Figure 9 - Unit 1 Measured Flow 



r i 


Figure 10 - Unit 1 Pressure and Temperature Interaction 


Noise at the low and high ends of the flow 
range was seen. The low end noise was caused by poor 
signal to noise ratios and electrical noise on very low 
signals. The high end noise was caused by flow 
turbulence and the inability to create a single flow 
passage capable of maintaining laminar flow over such 
a wide flow and pressure range. 

The concatenated error plot (Figure 11) shows 
that the concept of using two different flow 
measurement methodologies to cover the high and low 
flow ranges with a merging of the two methods in the 
middle was not capable of producing low errors in the 
merged area. This possibly could be corrected by use of 
multiple lookup tables instead of curve fitted equations 
but a superior approach appears to be that which is 
embodied in the by-pass unit which, uses only a single 
measurement methodology and avoids the merging of 
two separate range signals completely. Prior to actually 
testing the bypass unit a measure of its capability was 
obtained by studying the performance of unit 1 over a 
much narrower flow range (0.33 to 8 pounds per hour). 
Ignoring the large standard deviation due to noise, 
Figure 11 shows that the average error is almost 
^istant and capable of falling within a ±3% of reading 

£Bd. 


NITROGEN @ 100 PSIG AND 180 DEGREES F 



ACTUAL FLOW (Ib/hr) 

Figure 11 - Unit 1 Percentage Error 

Because of the averaging method used to 
reduce the noise levels the time delay associated with 
changing signals made it very difficult to distinguish 
between true hysteresis and time delay effects. .The 


8 



Figure 12 - Unit 1 Pressure and Y Attitude Interaction 


faster a sweep is run the wider the apparent hysteresis 
band due to the time delay. The best estimate of 
hysteresis is obtained by assuming that the narrowest 
band is representative of the hysteresis and is least 
affected by the time delay. A best fit straight line for 
increasing flow and another for decreasing flow were 
computed for the linear segment of the flow curve 
shown in figure 9. The hysteresis is half the difference 
between the two best fits and is approximately 0.7% of 
reading. The flow standard used has hysteresis of 0.1% 
of reading; thus, the unit 1 hysteresis is estimated to be 
0.6% of reading. Figure 13 shows a step response test 
for the sensor for one particular gas type, gas 
temperature, pressure and flow rate and is again 
representative of all of the data collected at other 
combinations of conditions. The sensor does not exhibit 
a typical linear system exponential time response, 
instead having a pure time delay of about 2.5 seconds 
before rising (or falling) instantly to the final value. 

The sensor uses a digital sampling and 

> raging technique to minimize the effect of flow noise 
I calculate the output mass flow signal, and thus 
nges the output once every 2 to 3 seconds. This low 



TIME (sec) 

UNIT 1 ACTUAL FLOW j 


Figure 13 - Unit 1 Step Response 

update rate is related to the use of non-optimized 
electronics in the unit and would be expected to be 


9 





increased in a production unit by the use of a more 
powerful micro-processor. 


Conclusions 



The existing unit can only be calibrated for a 
lie gas and has a flow range of 0.15 to 5 Ib/hr with an 
'accuracy of ± 0.15 Ib/hr. The goal for the new sensor 
was for a multi-gas capability and a range of 0.01 to 150 
Ib/hr with an accuracy of ± 3% of reading from 0.33 to 
150 Ib/hr and ± 0.01 Ib/hr from 0.01 to 0.33 Ib/hr. This 
goal represented a range multiplication of over 30 times 
the existing sensor and at 1 Ib/hr an accuracy 
improvement of a factor of 5. These were aggressive 
objectives. 

The test results for the type 1 sensor highlighted 
some specific problems with the unit but did confirm that 
the measurement principle is capable of operating over 
the flow, pressure and temperature ranges required. 
Furthermore, the testing verified that the unit is not 
affected by ambient pressure, temperature and humidity 
and by X-axis orientation. The problems identified by the 
testing were related to noise, accuracy, Y-axis 
orientation and unit to unit manufacturing repeatability. 
Preliminary evaluation of the data indicates that the 
problems with noise and accuracy can be alleviated and 
the problem with Y-axis orientation is not believed to be 
severe. However, the problem with unit to unit variation 
severely reduces the value of the unit by preventing 
interchangability. The problem related to accuracy is 
largely related to problems in the mid-flow range 

• )ciated with the merging of signals generated by the 
different measurement methodologies. 

The testing with the type 2 sensor is not 
sufficiently advanced to be included in the paper at this 
time but it is intended to include the data in an update of 
the paper. Nonetheless, the preliminary testing 
indicates, as expected, that this unit provides better 
accuracy and less sensitivity to orientation than the type 
1 unit. The preliminary testing of the type 2 sensor which 
has been accomplished has covered a flow range of 

0.01 to 150 Ib/hr for oxygen, nitrogen and helium at 
pressures from 50 to 900 psig and temperatures from 30 
to 150 degrees Fahrenheit. Generally from 2 Ib/hr up to 
150 Ib/hr the unit accuracy is very good. 

The test program has shown that the type 1 
sensor probably will not be able to completely reach the 
goals which were set. But with correction of a 
temperature sensor isolation problem and improved 
calibration at the low flow end of the range the type 2 
unit may be able to achieve the target objectives. 
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UNIT # 1 NITROGEN FLOW AT 70 DEGREES F 



(JM/qi) mou iiNn 


ACTUAL FLOW (Ib/hr) 



UNIT # 1 NITROGEN AT 70 DEGREES F 



ACTUAL FLOW (Ib/hr) 



UNIT 1 


DOE 1 


RUN 

NUMBER 

FLOW 

ACTUAL 

FLOW 

RS232 


FLOW 

ANALOG 

ERROR 
% READING 







1 

5.06 

5.307 

4.88 

5.356 

5.85 

2 

50.6 

54.2 

7.11 

53.3 

5.34 

3 

5.02 

5.375 

7.07 

5.39 

7.37 

BB 

50.65 

54.2 


52.89 

4.42 | 

BH 

5.015 

5.279 

■ 

5.397 


SB 

50.14 

57.53 

KgS 




5.074 

5.245 



— ' w 

2 

49.36 

65.7 

33.1 


32.66 | 


Table 1 
































File : 
Title: 
Cr—”nent : 

r^Bfcpnse : 




FL0W1 . DAT 

Solid State Flow Sensor Screening Test 
Screening to Reduce Number of Factors 
Flow & Err, Pwr, Delta P & Thermal Err 


Ax— j.ysis Of Variance (ANOVA) Table 


# 

Factor 

DOF 

SOS 

Percent 

Mean-Sqr 

Fcalc 

% Conf 

Effect 

1 

Flow Rate 

1 

5534 . 78 

98.43 

5534.783 

647.4 

97.50 

52.6060 


Ln 

1 

5534.78 

98.43 


647.4 



2 

Pressure 

1 

8.41357 

0.15 

8.413574 

. 9841 

-1.000 

2.05100 


Ln . 

1 

8.41318 

0.15 


. 9841 

- 


3 

Gas Temp 

1 

26.9097 

0.48 

26.90967 

3.148 

67.32 

3 .66800 


Ln 

1 

26.9084 

0.48 


3 . 147 



4 

Air Press 

1 

8.27490 

0.15 

8.274902 

.9679 

-1.000 

2 . 03400 


Ln 

1 

8.27430 

0.15 


.9678 



5 

Humidity 

1 

28.0806 

0.50 

28 . 08057 

3.285 

67.90 

3.74700 


Ln 

1 

28.0800 

0.50 


3.284 



6 

X- Attitude 

1 

8.13721 

0.14 

8.137207 

.9518 

-l.OOOi 

2.01700 


Ln 

1 

8.13656 

0.14 


.9517 



M 

Residual 

1 

8.54932 

0.15 

8.549316 




R 

Replicate 

0 

.000000 

0.00 

.0000000 





Pooled 

1 

8.54932 

0 . 15 

8.549316 





Total 

7 

5623.15 

100 . 00 





SOS 

Grand Total : 


13613.9 






SOS 

Mean (CF) : 


7990 . 75 






SOS 

Total 


5623.15 







Table 2 


File : 
Title : 
Comment : 
rations e 


FL0W1 . DAT 

Solid State Flow Sensor Screening Test 
Screening to Reduce Number of Factors 
Flow & Err, Pwr, Delta P & Thermal Err 


A^^ysis Of Variance (ANOVA) Table 


# 

Factor 

DOF 

SOS 

Percent 

Mean-Sqr 

Fcalc 

Fcrit 

Effect 

1 

Flow Rate 

1 

214.038 

31.77 

214.0382 

3.370 

161.4 

10.3450 


Ln 

1 

214.038 

31.77 

- 

3.370 



2 

Pressure 

1 

43.0593 

6.39 

43.05933 

.6780 

161.4 

4.64000 


Ln ' 

1 

43.0592 

6.39 


.6780 


- 

3 

Gas Temp 

1 

115.520 

17.14 

115.5200 

1.819 

161.4 

7 .6000 0 


Ln 

1 

115.520 

17.14 


1.819 



4 

Air Press 

1 

40.3203 

5.98 

40.32031 

.6349 

161.4 

4.49000 


Ln 

1 

40.3202 

5.98 


.6349 



5 

Humidity 

1 

171.495 

25.45 

171.4954 

2.700 

161.4 

9.26000 


Ln 

1 

171.495 

25.45 


2.700 



6 

X- Attitude 

1 

25.8482 

3.84 

25.84821 

.4070 

161.4: 

3.59500 


Ln 

1 

25.8480 

3 . 84 


.4070 



M 

^Residual 

1 

63.5058 

9.43 

63 .50580 





Pleplicate 

0 

.000000 

0.00 

.0000000 





Pooled 

1 

63.5058 

9.43 

63.50580 





Total 

7 

673 . 787 

100.00 






SOS Grand Total: 1525.39 
SOS Mean (CF) : 851.606 
SOS Total : 673.787 


Table 3 



UNIT 1 


D0E1A 


RUN 

NUMBER 

FLOW 

ACTUAL 

FLOW 

RS232 

ERROR 
% READING 

FLOW 

ANALOG 

ERROR 
% READING 







1 

5.02 

5.355 

6.67 

5.36 

6.77 

2 

50.79 

54.069 

6.46 

53.74 

5.81 

3 

5.0556 

5.438 

7.56 

5.429 

7.39 

4 

50.21 

52.05 

3.66 

51.66 

2.89 


Table 4 






























File : 
Title: 


Crrment : 

1 ^“vonse : 


FLOW1A .DAT 

Solid State Flow Sensor Screening Test 
Screening to verify humidity effect 
Flow & Err 


D^^Tgn Listing: 


Design Name: 
Design: 

# Treatments: 


[L4 ] 

Full Factorial 
C 4 ] 


[2"k] , (K = 2 ) 


Extra' DOF : 


[ 0] , Extra DOF For > 2 -Factor Ints 


TC# 


1 2 


1: 
2 : * 
3 : 



FLOW 

HUMIDITY 

5 

25 

50 

25 

5 

85 

50 

85 

FLOW 

HUMIDITY 


Conf ounding/Aliasing Information : 


No Confounding Exists Since P=0 


Table 5 



File : 
Title : 
Comment : 
F^fcpnse : 




FL0W1A . DAT 

Solid State Flow Sensor Screening Test 
Screening to verify humidity effect 
Flow. Sc Err £S2_32L 


A*.^ysis Of Variance (ANOVA) Table 


# Factor 

DOF 

SOS 

Percent 

Mean-Sqr 

Fcalc 

Fcrit 

Effect 

1 FLOW 

1 

2237.81 

99.90 

2237.811 

1939. 

161.4 

47.3055 

Ln 

1 

2237.81 

99.90 


1939. 



2 HUMIDITY 

1 

1.01147 

0.05 

1.011475 

.8765 

161.4 

-1.00550 

Ln. 

1 

1.01103 

0.05 


.8761 

- 


Residual 

1 

1.15405 

0.05 

1.154053 




Replicate 

0 

.000000 

0.00 

. 0000000 




Pooled 

1 

1.15405 

0.05 

1.154053 




Total 

3 

2239.98 

100.00 






SOS Grand Total: 5614.95 
SOS Mean (CF) : 3374.97 
SOS Total : 2239.98 



Table 6 


File : 
Title: 
Cement : 
>nse : 




FL0W1A . DAT 

Solid State Flow Sensor Screening Test 
Screening to verify humidity effect 
Flow & Err f?S2.12_ 


Al^^ysis Of Variance (ANOVA) Table 


Factor 

DOF 

SOS 

Percent 

Mean-Sqr 

Fcalc 

% Conf Effect 

FLOW 

1 

4.22301 

49.46 

4.223007 

1.241 

53.42 -2.05500 

Ln 

1 

4.22302 

49.46 

' 

1.241 


HUMIDITY 

1 

.912018 

10.68 

.9120178 

.2679 

-1.000 -.955000 

Ln . 

1 

.912025 

10.68 


.2679 

- 

Residual 

1 

3.40404 

39.86 

3.404037 


* • 

Replicate 

0 

.000000 

0.00 

.0000000 



Pooled 

1 

3.40404 

39.86 

3.404037 



Total 

3 

8.53906 

100.00 





SOS Grand Total: 156.770 
SOS Mean (CF) : 148.231 
SOS Total : 8.53906 



Table 7 

















r i 

L A 


f ^ 

L A 


f ^ 

L A 



Table 8 










































File : 
Title: 
Ov'Tient : 
mse 




FLOW_LOG . DAT 

Solid State Flow Sensor Screening Test 
Preliminary Sensitivity Identification 
Flow Sc Err, Pwr Led Cob£b 


An~xysis Of Variance (ANOVA) Table 


# 

Factor 

DOF 

SOS 

Percent 

Mean-Sqr 

Fcalc 

Fcrit 

1 

Gas 

2 

. 003643 

0 . 03 

. 0018215 

.3176 

4.459 


Ln 

1 

.001976 

0 . 02 

- 

.3446 



Q'd 

1 

. 001676 

0.01 


.2922 


2 

Flow 

2 

12.5396 

99.43 

6.269806 

1093 . 

4.459 


Ln 

1 

12.5371 

99.41 


2186. 



Qd 

1 

. 002544 

0.02 


.4435 


3 

Press 

2 

.021526 

0.17 

.0107632 

1.877 

4.459 


Ln 

1 

.006238 

0 . 05 


1.088 



Qd 

1 

. 015293 

0.12 


2.667 


4 

GasTemp 

2 

.000948 

0 . 01 

. 0004740 

. 0826 

4.459 


Ln 

1 

. 000468 

0 . 00 


. 0815 



Qd 

1 

. 000488 

0 . 00 


. 0850 


5 

Y Attitude 

1 

.000128 

0 . 00 

.0001278 

. 0223 

5.318 

m 

| Ln 

1 

.000134 

0 . 00 


. 0234 


m 

^Residual 

8 

.045879 

0.36 

. 0057349 




Replicate 

0 

. 000000 

0 . 00 

.0000000 




Pooled 

8 

. 045879 

0.36 

. 0057349 




Total 

17 

12.6117 

100 . 00 




SOS 

Grand Total : 


33 . 0845 





SOS 

Mean (CF) 


20.4728 





SOS 

Total : 


12.6117 






Signif . 


Table 9 


File: FLW_LG_E . DAT 

Title: Solid State Flow Sensor Screening Test 
Co-mient: Preliminary Sensitivity Identification 
F^fepnse: Flow & Err , Pwr £s 2 . 5 z Hector 




A.. .^ysis Of Variance (ANOVA) Table 


# 

Factor 

DOF 

SOS 

Percent 

Mean-Sqr 

Fcalc 

Fcrit 

1 

Gas 

2 

31.7146 

4.19 

15.85730 

.3015 

4.459 


Ln 

1 

24.4245 

3.23 

- 

.4643 



Qd 

1 

7.29000 

0.96 


.1386 


2 

Flow 

2 

99.2856 

13.12 

49.64282 

.9437 

4.459 


Ln 

1 

90.2557 

11.92 


1.716 



Qd 

1 

9.03002 

1.19 


.1717 


3 

Press 

2 

196.264 

25.93 

98.13184 

1.866 

4.459 


Ln 

1 

48.8840 

6.46 


.9293 



Qd. 

1 

147.380 

19.47 


2 . 802 


4 

GasTemp 

2 

6.60510 

0.87 

3.302551 

.0628 

4.459 


Ln 

1 

.330008 

0.04 


.0063 



Qd 

1 

6.27502 

0.83 


.1193 


5 

Y Attitude 

1 

2.30414 

0.30 

2.304138 

.0438 

5.318 

m 

^Ln 

1 

2.30409 

0.30 


.0438 


% 

Residual 

8 

420 . 820 

55.59 

52.60254 




Replicate 

0 

.000000 

0.00 

.0000000 




Pooled 

8 

420.820 

55.59 

52.60254 




Total 

17 

756.993 

100.00 




SOS 

Grand Total : 


1421.66 





SOS 

Mean (CF) : 


664.666 





SOS 

Total : 


756.993 






Signif . 


Table 10 















Flow Sensor Screening T 

iary Sensitivity Identificoti 

>r, Pwr Ljo& &S2-V 

























UNIT 1 DOE 3A 


RUN 

NUMBER 

GAS 

TEMP 

GAS 

PRESS 

FLOW 

ACTUAL 

FLOW 

RS232 

CORR 

RS232 

ERROR 
% READING 

FLOW 

ANALOG 

CORR 

ANALOG 

ERROR 
% READING 



1SBBB1|| 





■MM8S 

■■ 


1 

40 

70 

1.99 

2.182 

2.1930 

9.65 

2.195 

2.2060 

10.30 

2 

165 

70 

1.989 

2.172 

2.1840 

9.20 

2.163 

2.1750 

8.75 

3 

40 

850 

2.024 

2.251 

2.2243 

WBB3M 

2.265 

2.2381 

11.91 

4 

165 

850 

1.9928 

2.172 

2.1798 

8.99 

2.165 

2.1728 

8.64 


Table 11 


UNIT 1 DOE 3B 


RUN 

NUMBER 

GAS 

TEMP 

Y 

ATTDE 

FLOW 

ACTUAL 

FLOW 

RS232 

CORR 

RS232 


FLOW 

ANALOG 

CORR 

ANALOG 



iiiiiSifl 







wmmi 


1 

40 

UP 

2.014 

2.266 

2.2502 

12.51 

2.282 

2.2661 

13.31 

.2 

165 

UP 

1.991 

2.233 

2.2431 

12.15 

2.234 

2.2441 

12.20 

3 

40 

DOWN 

2.006 

2.152 

2.1456 

7.28 

2.163 

2.1565 

7.83 

4 

165 

DOWN 

1.986 

2.152 

2.1672 

8.36 

2.155 

2.1702 

8.51 


Table 12 

UNIT 1 DOE 3C 


RUN 

NUMBER 

GAS 

PRESS 

Y 

ATTDE 

FLOW 

ACTUAL 

FLOW 

RS232 

CORR 

RS232 

ERROR 
% READING 

FLOW 

ANALOG 

CORR 

ANALOG 

ERROR 
% READING 











1 

100 

UP 

2.0166 

2.222 

2.2037 

10.19 

2.239 

2.2206 

11.03 

2 

850 

UP 

2.014 

2.266 

2.2502 

12.51 

2.282 

2.2661 

13.31 

3 

100 

DOWN 

2.008 

2.207 

2.1982 

9.91 

2.221 

2.2122 

10.61 

4 

850 

DOWN 

2.006 

2.152 

2.1456 

7.28 

2.163 

2.1565 

7.83 


Table 13 


UNIT 1 DOE 3D 


RUN 

NUMBER 

GAS 

PRESS 

GAS 

TYPE 

FLOW 

ACTUAL 

FLOW 

RS232 

CORR 

RS232 

ERROR 
% READING 

FLOW 

ANALOG 

CORR 

ANALOG 

ERROR 
% READING 











1 

70 

N2 

1.989 

2.172 

2.1840 

9.20 

2.163 

2.1750 

8.75 

2 

850 

N2 

1.9928 

2.172 

2.1798 

8.99 

2.165 

2.1728 

8.64 

3 

70 

He 

2.098 

3.064 

2.9209 

46.04 

3.067 

2.9237 

46.19 

4 

850 

He 

1.996 

2.1767 

2.1811 

9.05 

2.177 

2.1814 

9.07 


Table 14 


UNIT 1 DOE 3E 


RUN 

NUMBER 

GAS 

TEMP 

GAS 

TYPE 

FLOW 

ACTUAL 

FLOW 

RS232 

CORR 

RS232 

ERROR 
% READING 

FLOW 

ANALOG 

CORR 

ANALOG 

ERROR 
% READING 











1 

40 

N2 "1 

1.99 

2.182 

2.1930 

9.65 

2.195 

2.2060 

10.30 

2 

165 

N2 

1.989 

2.172 

2.1840 

9.20 

2.163 

2.1750 

8.75 

1 3 

40 

He 

1.9917 

2.9597 

2.9720 

48.60 

2.975 

2.9874 

49.37 

' 4 

165 

He 

2.098 

3.064 

2.9209 

46.04 

3.067 

2.9237 

46.19 


Table 15 







































































































































































































































File : 
Title : 
C^-~?ent : 
lse : 




FL0W3_A .DAT 

Solid State Flow Sensor Screening Test 
Identification of Interactions 
Flow & Err 


An-xysis Of Variance (ANOVA) Table 


# 


1 


2 


Factor 

DOF 

SOS 

Percent 

GAS TEMP 

1 

.000719 

59.09 

Ln 

1 

.000716 

58.80 

GAS PRESS 

1 

.000185 

15.20 

Ln. 

1 

.000184 

15.09 

Residual 

1 

.000313 

25.71 

Replicate 

0 

.000000 

0 . 00 

Pooled 

1 

. 000313 

25.71 

Total 

3 

.001217 

100.00 


Mean-Sqr 

Fcalc 

Fcrit 

Effect 

.0007191 

2.299 

2.288 

161.4 

- .026750 

.0001850 * 

.5915 

.5870 

161.4 

.013550 

. 0003128 
. 0000000 
. 0003128 





SOS Grand Total: 19.2781 
SOS Mean (CF) : 19.2769 
SOS Total : .001217 



Table 16 
















File: FL0W3_B .DAT 

Title: Solid State Flow Sensor Screening Test 

-nent : Identification of Interactions 

I ^^nse : Flow & Err 


i m 

AiiSWy, 


sis Of Variance (ANOVA) Table 


# 


1 


2 


Factor 

DOF 

SOS 

Percent 

GAS TEMP 

1 

.000055 

0.66 

Ln 

1 

.000053 

0.63 

Y Attitude 

1 

.008148 

96.96 

Ln . 

1 

.008145 

96.92 

Residual 

1 

.000200 

2.38 

Replicate 

0 

.000000 

0.00 

Pooled 

1 

.000200 

2.38 

Total 

3 

.008404 

100.00 


Mean-Sqr 

Fcalc 

Fcrit 

Effect 

.0000553 

.2762 

.2625 

161.4 

.007250 

.0081482 

40.69 

40.67 

161. .4 

-.090250 

.0002003 

.0000000 

.0002003 





SOS 

Grand Total : 

19.3953 

SOS 

Mean (CF.) : 

19.3868 

SOS 

Total : 

.008404 



Table 17 





















File : FL0W3 

C .DAT 






Title: Solid 

State 

Flow Sensor Screening Test 




O — -sent: Identification of Interactions 


- 



: g 

^ise: Flow 

& Err 







J 

^sis Of Variance 

(ANOVA) Table 





# 

Factor 

DOF 

SOS 

Percent 

Mean-Sqr 

Fcalc 

Fcrit 

Effect 

1 

GAS PRESS 

1 

.000006 

0.10 

.0000057 

.0023 

161.4 

-.003050 


Ln 

1 

.000009 

0.17 


.0038 



2 

Y Attitude 

1 

.003027 

55.12 

.0030270 

1.231 

161.4 

-.055050 


Ln. 

1 

.003031 

55.19 


1.233 


- 


Residual 

1 

.002459 

44.77 

.0024586 





Replicate 

0 

.000000 

0.00 

.0000000 





Pooled 

1 

.002459 

44.77 

.0024586 





Total 

3 

.005491 

100.00 





SOS 

Grand Total 


19.3554 






SOS 

Mean (CF-V 


19.3499 





i 

SOS 

Total 


.005491 








Table 18 














File : 
Title : 

C • T nent : 
inse : 


AiSffy 


FLOW3JO .DAT 

Solid State Flow Sensor Screening Test 
Identification of Interactions 
Flow Sc Err 


Al^ffysis Of Variance (ANOVA) Table 


# 


1 


2 


Fact ox 

DOF 

SOS 

Percent 

GAS PRESS 

1 

.138384 

33.76 

Ln 

1 

.138384 

33.76 

GAS TYPE 

1 

.136236 

33.24 

Ln. 

1 

.136235 

33.24 

Residual 

1 

.135275 

33.00 

Replicate 

0 

.000000 

0.00 

Pooled 

1 

.135275 

33.00 

Total 

3 

.409895 

100.00 


Mean-Sqr 

Fcalc 

Fcrit 

Effect 

.1383839 

1.023 

1.023 

161.4 

- .372000 

.1362362 

1.007 

1.007 

161.4 

.369100 

.1352749 

.0000000 

.1352749 





SOS Grand Total: 22.8102 
SOS Mean (CE) : 22.4003 
SOS Total : .409895 



Table 19 




















File : 
Title: 
Comment : 


a 


FL0W3_E .DAT 

Solid State Flow Sensor Screening Test 
Identification of Interactions 
Flow St Err 


Ai^ffysis Of Variance (ANOVA) Table 


# 


1 


2 


Factor 

DOF 

SOS 

Percent 

GAS TEMP 

1 

.000902 

0.16 

Ln 

1 

.000903 

0.16 

GAS TYPE 

1 

.574484 

99.77 

Ln . 

1 

.574488 

99.77 

Residual 

1 

.000446 

0 . 08 

Replicate 

0 

. 000000 

0.00 

Pooled 

1 

.000446 

0.08 

Total 

3 

.575832 

100.00 


Mean-Sqr 

Fcalc 

Fcrit 

Effect 

.0009022 

2.021 

2.023 

161.4 

- .030050 

.5744839 

1287. 

1287. 

161.4 

.757950 

. 0004463 
. 0000000 
.0004463 





SOS Grand Total: 26.9435 
SOS Mean (CF) : 26.3677 
SOS Total ‘ : .575832 


Table 20 


















UNIT 1 DOE 3D REPEAT 


RUN 

NUMBER 

GAS 

PRESS 

GAS 

TYPE 

FLOW 

ACTUAL 

FLOW 

RS232 

CORR 

RS232 

ERROR 
% READING 

FLOW 

ANALOG 

CORR 

ANALOG 

ERROR 
% READING 

* „ 










1 

100 

N2 

5.017 

5.607 

5.5880 

11.76 

5.639 

5.6199 

■KE2Z3HI 

2 

850 

N2 

5.051 


5.5217 

10.43 

5.598 

5.5415 

■BilXI 

3 

100 

He 

0.5027 

0.5298 

5.2695 

5.39 

0.5317 

5.2884 

5.77 

4 

850 

He 

0.4978 

0.5145 

5.1677 

3.35 

0.5154 

5.1768 

3.54 


Table 21 
















































File: FL0W3_D .DAT 

Title: Solid State Flow Sensor Screening Test 

Cement: Identification of Interactions 

onse : Flow & Err 




sis Of Variance (ANOVA) Table 


# Factor 

DOF 

SOS 

Percent 

Mean-Sqr 

Fcalc 

Fcrit 

Effect 

1 GAS PRESS 

1 

.007065 

5.87 

.0070648 

21.53 

161.4 

-.084050 

Ln 

1 

.007064 

5.87 


21.53 



2 GAS TYPE 

1 

.113052 

93.86 

.1130524 

344.6 

161.4 

-.336250 

Ln. 

1 

.113064 

93.87 


344.6 

- 


Residual 

1 

.000328 

0.27 

.0003281 




Replicate 

0 

.000000 

0 . 00 

. 0000000 




Pooled 

1 

.000328 

0.27 

.0003281 




Total 

3 

.120445 

100 . 00 






SOS Grand Total: 116.188 
SOS Mean (CF t ) : 116.067 
SOS Total : .120445 



Table 22 











Title. Solid State Flow Sensor Screening Test ' File:: FLOW3_D.DAT 

Comment: Identification of Interactions Date: 03-1 5- 1 996 











HONEYWELL FLOW SENSOR UNIT 1 


GAS 

PRESS 


ERROR (0.33 TO 8 LB/HR) | 

AVERAGE 

STD DEV 

MAX 

MIN 









900 

30 

8.6 

7.15 


-14.2 

He 

100 

30 

8.1 

6.61 

76.2 

-5.9 

He 

900 

180 

8 

4.94 

40.1 

-7.3 

He 

100 

180 

8.8 

5.68 

74 

-9.7 

N2 


30 

11.1 

HHEES 

40.1 

-8 

N2 

100 

30 

10.5 

■HES 

59.8 

-0.5 

N2 


180 

9.4 

5.3 

39.1 

-16.8 

N2 

100 

180 

10.9 

3.41 

34.6 

-2.4 

02 

900 

30 

13.1 

3.3 


1.5 

02 

100 

30 

9.5 

5.06 

36.3 

-24.2 

02 

900 

180 

8.3 

4.78 

32.8 

-17.4 

02 

100 

180 

9.8 

4.67 

25.3 

-15.5 


Table 23 
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HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 







HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 








HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 




HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (IbAir) 









HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @100 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 
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HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ibfttr) 



HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 



HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 900 PSIG AND 30 DEGREES F) 



Umi) Mold aaunsvaw tihmabnoh 


ACTUAL FLOW (Ib/hr) 
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HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 




HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 
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HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 





HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ih/hr) 









HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (lb/hr) 



HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 180 DEGREES F) 



(JM/qi) MOU Q3unsvaw TI3MA3NOH 


ACTUAL FLOW (Ib/hr) 
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HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 


HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 180 DEGREES F) 



(jq/qi) mou aaunsvaw ih3ma3noh 


ACTUAL FLOW (Ib/hr) 







HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 










HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 






HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 




HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 180 DEGREES F) 



(JM/qi) MOU Q3UnSV3W 113MA3NOH 


ACTUAL FLOW (Ib/hr) 



HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 










HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 100 PSIG AND 180 DEGREES F) 



(•iq/qi) MOid aaunsv3w thmasnoh 


ACTUAL FLOW (Ib/hr) 
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HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 900 PSIG AND 180 DEGREES F) 



On/qi) mou aaansv3w TI3MA3NOH 


O 


ACTUAL FLOW (lb/hr) 










HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 900 PSIG AND 180 DEGREES F) 



(JM/Ql) MOId 03UnSV3W TI3MA3NOH 


O 












HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 900 PSIG AND 180 DEGREES F) 



(JM/QI) MOId 03UnSV3W USMASNOH 


O 


ACTUAL FLOW (Ib/hr) 













HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(NITROGEN @ 900 PSIG AND 180 DEGREES F) 



OM/qi) Mou aaynsvaiN usmasnoh 


O 


ACTUAL FLOW (Ib/hr) 
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HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 100 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 
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HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 100 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 











HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 100 PSIG AND 30 DEGREES F) 



(jn/qi) mou aaunsvaw naMA3NOH 


ACTUAL FLOW (Ib/hr) 
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HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 





HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 30 DEGREES F) 


O 



ACTUAL FLOW (Ib/hr) 




HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 30 DEGREES F) 


o 



ACTUAL FLOW (Ib/hr) 







HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 






HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSSG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 





HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 



HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 30 DEGREES F) 



(JM/qi) MOU aaunSV3W T13MA3NOH 


ACTUAL FLOW (Ib/hr) 
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HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 







HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 




HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 








HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 100 PSIG AND 180 DEGREES F) 


CO 



ACTUAL FLOW (lb/hr) 




















HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 


HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 
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HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 180 DEGREES F) 



(jn/qi) Mend aaunsvaw na/WA3NOH 


ACTUAL FLOW (Ib/hr) 



HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (tb/hr) 



HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 180 DEGREES F) 



L A 


HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (lb/hr) 



90*0 
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HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 180 DEGREES F) 


CO 



ACTUAL FLOW (Ib/hr) 














HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (lb/hr) 




HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 











HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(OXYGEN @ 900 PSIG AND 180 DEGREES F) 


0 



ACTUAL FLOW (lb/hr) 








rrri 




HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 100 PSIG AND 30 DEGREES F) 



Figure 106 







(6u|peaj jo %) U0UU3 MOId 113MA3N0H 



HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 100 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (lb/hr) 












HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 100 PSIG AND 30 DEGREES F) 


0 



ACTUAL FLOW (Ib/hr) 






HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @100 PSIG AND 30 DEGREES F) 



Figure 1 1 1 





HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 100 PSIG AND 30 DEGREES F) 



Figure 112 







HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @100 PSIG AND 30 DEGREES F) 



Figure 113 






HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 100 PSIG AND 30 DEGREES F) 



Figure 114 







HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 


HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 




HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 








HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 






HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 





HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 



HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 900 PSIG AND 30 DEGREES F) 
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ACTUAL FLOW (Ib/hr) 





HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 












HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 900 PSIG AND 30 DEGREES F) 



ACTUAL FLOW (Ib/hr) 









HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 



HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 



HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 100 PSIG AND 180 DEGREES F) 



(Bu|pe»jjo%) UOUU3 MOU T13MA3NOH 


ACTUAL FLOW (lb Air) 


HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 100 PSIG AND 180 DEGREES F) 



(6u|peaj jo %) UOUU3 MOU TI3MA3NOH 


ACTUAL FLOW (Ib/hr) 
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HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 100 PSIG AND 180 DEGREES F) 



(jq/qi) Mold aaunsvam tismasnoh 



HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 100 PSIG AND 180 DEGREES F) 



ACTUAL FLOW (Ib/hr) 







HONEYWELL UNIT 1 - MEASURED FLOW VERSUS ACTUAL FLOW 
(HELIUM @ 100 PSIG AND 180 DEGREES F) 
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GAS FLOW SENSORS FOR SPACE APPLICATIONS 
Phase I. Fabrication of Two Demonstration Units 

1 . EXECUTIVE SUMMARY 

The purpose of this program was to demonstrate the performance of a 
thermal flow microsensor. It would be intended for the measurement of 
flows of helium, nitrogen, oxygen and hydrogen in both the ARPCS 
(Atmospheric Revitalization Pressure Control Subsystem) and in the FCS 
(Fuel Cell System) , which have been using sensors that experienced high 
failure rates. 

An existing sensor design was adapted for this project. Two sensors were 
built and tested, and delivered to Car let on for further tests, after 
improvements were made to overcome interferences by what was recognized 
as thermal microconvection. 

In spite of -the improvements, the test results indicate that not all 
performance attributes materialized, which had been predicted for *■ 
operation at pressures up to 61 bar (900 psi) , based on previous data of 
flow microsensor performance up to 4 bar (60 psia) . The unforeseen 
factor that increasingly interfered as pressure was increased, was one 
caused by thermal microconvection currents generated by the sensor's 
microheaters. These caused both large, signal shifts as well as 
deterioration of the S/N ratio. Whereas the direct influence of these 
currents could be reduced significantly via implemented improvements < 14 ) , 
the ability of a single point flow sensor to meet the accuracy 

• specifications over the large 15,000:1 dynamic flow range remained 
marginal . 

Test results of these improved units, while meeting expectations in 
terms of the response time, power consumption, insensitivity to 
orientation up to 500 psi, and accuracy at flows below 5 lb/h, were 
disappointing as far as accuracies achieved with the high flow sensing 
approach (above 5 lb/h) , labour-intensive adjustments of the matching of 
the high and low flow signals, flow signal errors resulting from the 
temperature difference between the flowing gas and the sensor block, in 
which the property and temperature sensors are mounted, the orientation 
sensitivity at 900 psig and signal noise, and the fact that the 
remaining orientation sensitivity also prevented us from differentiating 
between N 2 and O 2 . 

We therefore recommended that the direct, single-point flow measurement 
of Units#l and 2 be replaced by a bypass, multiple-point flow sensing 
approach. Such recommendation was made on the basis of tests made with 
two such units, which not only showed predictable gas composition 
correction behaviour, a possibility to replace the heater power approach 
used at high flows by the more accurate differential approach over the 
the full flow range, but also a more linear and predictable output, 
enhanced S/N and the possibility of further reducing the effect of 
orientation-dependent, thermal microconvection, due to the beneficial 
influence of mounting the sensor chip in a well-defined, more 
^^constraining microchannel . 
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2 . INTRODUCTION 


The purpose of this program was to adapt an existing design, build, 
calibrate and test two sensors for the measurement of gas mass flow of 
helium, nitrogen, oxygen and hydrogen; this single design flow sensor 
was to serve in multiple applications in both the ARPCS (Atmospheric 
Revitalization Pressure Control Subsystem) and in the FCS (Fuel Cell 
System) , which have been using sensors that experienced high failure 
rates. The redesigned sensors were to be able to meet the target 
specifications listed in Table 1, with the noted exceptions for the 
demonstration units completed in this phase of the program. Ultimately 
the redesigned production sensors are to replace those shown in Figs. la 
and b. 


The main deliverables were to consist of: 

- Two testable flow sensors, with all fluid and electrical mating 
connectors needed to interface with existing NASA equipment, with* 
circuitry* to provide a 0-5 V output representing a compensated flow 
range from 0.1 to 150 lb/h; and 

- This final report, documenting all applicable performance and accuracy 
evaluation test results, in relation to the requirements listed in 
Table 1. 


The Honeywell effort was to include flow testing with He (as a simulant 
of H 2 ) , O 2 and N 2 at pressures between 200-900 psia, up to mass flows of 
15 lb/h, and consulting support on tests at Carleton and at NASA- JSC. 


3. TECHNICAL DISCUSSION 




1. Original Design and Fabrication Plan: Honeywell-HTC was to 
adapt its compensated flow microsensor design, see Fig. 3, to NASA- 
JSC's application. As agreed with Carleton, both deliverable units 
were to use the connectors and fittings shown in Figs.l and 2, to be 
provided by Carleton. The main features of this plan involved tasks 
to: 

- Use existing gas flow microsensor chips 

- Adapt existing analog front-end and digital linearization and 
compensation circuit designs and printed circuit boards; 

- Mate the inlet and outlet fittings to existing flow channel 
designs, which have been used in developmental, electronic gas 
meters to sense flows ranging from less than 0.01 to over 150 
lb/h with one single design; 

- Design a sealed housing for the electronics, with provisions to 
accept the specified connectors and sensor inputs; and 

- Use and/or adapt Honeywell's T, p and gas composition correction so 
that one single gas mass flow sensor design can be used with any of 


the required gases. 

- Calibrate each flow sensor system with N 2 , O 2 and He, at atmospheric 
pressure and covering the specified temperature range, with mass 
flows from 0.01 to 15 lb/h; and with air up to 150 lb/h. 

- Provide consulting help in support of comprehensive tests at ,■ 


Carleton (up to 900 psia) and at NASA with hydrogen. 
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A typical set of anticipated flow sensor performance data is shown in 
ig.4, obtained with a flow channel featuring a Venturi nozzle with a 14 
ID throat. For optimum performance in terms of mass flow range, low 
noise and output data accuracy we proposed to use the differential 
measurement approach at low flows, see e.g. Gd-u data in Fig. 4, and the 
heater-power approach for high flows, represented by the Htr. Power data 
in that same figure. At some intermediate range, e.g. 2-8 lb/h in Fig. 2, 
a weighted average between these two measurement approaches was to be 
formed for each sensor output. Fig. 4 also shows that the flow-induced 
pressure drop, even at atmospheric pressures (or about 60x below the 
maximum of 900 psia) , is about 3x below the performance goal of 28 psid, 
see Table 1, or llx further below that goal if the data had been taken 
at 200 psig. In view of that, and also to insure having (1) adequate 
signal level at low flows and (2) enough room to fit this flow channel, 
we agreed with the sponsors (Carleton and NASA) to reduce its cross 
section to half of its size, i.e. to reduce the throat diameter of the 
Venturi from 14 to 10 mm( 8) . As shown in Fig. 5, this reduction was not 
going to caiise excessive flow velocities even at pressures as low as, 40 
psia. 




At this meeting < 8 > we also reached tentative agreement about how to 
format the sensor's output in a way that would be compatible with the 
existing NASA software but not penalize the achievable accuracy of the 
new sensor. The agreed solution was to proportion the 0-5 volt output to 
a logarithmic, 0.001-150 lb/h mass flow range. This provides a 
graininess of only 1.94%, see Fig. 6, which could not be achieved with 
any of the linear scale options. Fig. 6a shows that this graininess is 
onsistent with the available bits (nine) and the desired dynamic range 
(15,000) . 


The 14 mm diameter sensor packages had been used before for: 

(1) developmental, natural gas metering applications W near atmospheric ( 
pressure, 

(2) up to 60 psia to prove that the mass flow sensor outputs were 
independent of pressure (2) , and 

(3) at about 3500 psia in a short field test with compressed natural 
gas; this verified that the sensor structure was strong enough to 
operate at such pressures, but also showed that further noise 
suppression work was needed. 

Having built the new, 10 mm ID Venturi, we checked that its performance 
was comparable to that of the previous 14 mm versions, see Fig. 7, by 
completing and plotting the data obtained at 20 and 50 psig. However, 
when the pressure was raised above 100 psig, unexpected signal 
distortions were observed, see Fig. 8, which became worse as the pressure 
was increased towards 900 psi. The output signal would drop in spite of 
maintaining the mass flow constant. The effects were large enough to 
dispel any hope of achieving a satisfactory sensor performance without 
major modifications. We therefore agreed with the sponsors that we 
needed to get back to fundamentals to fix this problem. Honeywell 
elected to fund this effort independently, because the solutions are 
needed also for other, non-space applications. 

j^.2 Thermal Flow Sensor Fundamentals: 

Standard volume or mass flow measurements are realized with two sensors 
that are positioned as indicated in Fig. 3: one measures flow while the 
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other senses thermal conductivity, specific heat, temperature and 
pressure. The latter is used to generate the compensation factor, C v , to 
obtain a fluid composition- and temperature-compensated standard volume 
flow, as detailed in refs. (1,11): 

V = V u /Cv, or m = m u /CM, or m = m u p/Cv 

with 

C V = (k/k Q ) nl (Cp/Cpo) n2 (Pr/Pr 0 ) n3 (T/T 0 ) n4 (1) 

where ni =-0.8827 or -0.8904 
n 2 = 1.1110 or 1.1500 

n 3 = 0 or 0.3003 

n 4 = 0.7000 or 0.6495< 13 > . 

The temperature exponent would ideally only depend on gas properties; in 
practice it' also is influenced by the temperature dependence of sensor 
materials and electronics. The use of one or the other set of exponents 
depends on the desired accuracy, and on whether on-line measurement of 
only k, c p and T, or k, c p , T| and T is available. As verified 
recently (11 ' 13> , these expressions for Cv are (1) universal for any type 
of thermal flow sensor, and (2) independent of pressure, within about 
±1%. It follows that for operation with known gases the Cv are known 
except for their temperature dependence, so that a system optimized for 
least cost and weight, and for max. reliability (i.e. min. number of 
parts) would not recompute all properties at every time interval, but 
rather: 

(a) receive an external input code to identify the gas being measured, 
so that the appropriate value of Cv can be activated from memory. 
Such a value, Cv=Cv 0 would be one corresponding to a standard 
temperature, T=T 0 , or 

(b) reduce the measurements to the minimum number required to identify 
the gas at hand, which would then be associated with the appropriate 
Cvo from memory. 

The sensor output signal for the low-flow region, AGi, can typically be 
represented by an equation close to V u i = a o +aiAGi 0 - 85 , with ao=0, see 
Fig. 9a, as long as the sensor signal, AGi=AGd-u/ represents a flow that 
corresponds to a temperature difference, AT^-u/ that is smaller than the 
heater temperature rise above the environmental temperature. At larger 
flows, the signal eventually saturates since ATd- u <AThtr • However, at 
such high flows, the heater power, G 2 , required to maintain the set AT h ti 
continues to rise, so that flows may be represented and determined by an 
expression, see also Figs. 4 and 7, of V u2 = b 0 +biG 2 ml +b 3 G 2 i:n2 , which is 
equivalent to the classical hot-wire anemometer equation, G 2 - H p , u2 - 
Hp,u2o • The above gas composition correction equations apply to it as 
well . 


The performance obtained with a flow channel of less than 70 mm (5»3 in.) 
in length, simply consisting of two screens and one honeycomb flow 
straightener upstream of a 14 mm ID Venturi nozzle, is illustrated in 
Fig. 9, where Fig. 9a shows the overall output signal shape of the 



uncompensated signal, for three gases, including N 2 , at flows below 15 
J-b/h . down to £5 g/h, i.e. £0.01 lb/h, and up to about 10 lb/h, before 
this differential measurement method saturates. Fig. 9b shows that the 
flow error after compensation for gas type and temperature is less than 
±3% of reading. 

In Fig. 4 the signal V<i_ u saturates at about 25 lb/h; however, flows 
greater than 15 lb/h and even beyond 150 lb/h can be measured with the 
same flow channel, without undue flow noise, by sensing the power 
required by the sensor heating element to maintain its temperature at a 
constant temperature rise above ambient (~100°C); while these 
measurements were taken at atmospheric pressure, no increase in flow 
noise is expected with flows at 1000 psia since the Reynolds number is 
pressure independent at constant mass flow. Fig. 4 also shows that the 
pressure drop at 150 lb/h air flow, (circle points) , was less than 3 
psid. 

Fig. 10 shows that gas composition or "type of ga s" compensation can well 
be provided within the ±3% error margin for gases differing in ' 

properties a's much as H 2 , He, N 2 , air and Ar-CH 4 mixtures; this set o.f 
measurements was obtained with a commercial AWM2100V type of Honeywell 
sensor, but could have been obtained with the above flow channel as 
well. 


The composition and temperature correction is performed by processing 
the outputs of the flow sensor and those of the property sensor, as 
described in ref . (1) , with an Intel 87C51 microcontroller; Fig.2e 
graphically illustrates how thermal conductivity and specific heat are 
pleasured. Flow measurements of up to 4 bar (60 psia) had indicated that 
pressure changes do not influence the flow sensor output, as would be 
xhe case for a true mass flow sensor. However, at pressures of up to 
1000 psia, real gas changes in compressibility factor, thermal 
conductivity and specific heat (which for ideal gases are pressure- 
independent) may require that pressure be measured and entered in a 
compensation scheme. Should that be necessary, the property sensor's 
built-in Pirany-type pressure sensor will be called into service, as 
described before (1 ' 7) ; its location on the chip is pointed out in the 
microphotographs of Fig.ll <7) . We have used this inherently burst-proof 
pressure sensor for pressure measurements of up to 2000 psia< 7 >. Thermal 
conductivity and specific heat measurements have been routinely carried 
out to 60 psig, and more recently to 600 psia in our automated 
calibration system. 


The weight s and volume s of the demonstration units were estimated by 
adding the contributions of the housing, which is assumed to be made of 
aluminum (density of 2.702 g/cm 3 ) for two single-sided circuit boards of 
75 x 150 mm and those of the flow channel of 3.6 mm OD. By subtracting 
the inner volumes, the weight of the demonstration unit is 0.75 lb, thus 
leaving 0.6 lb for electronics, fittings, connectors, sensors and flow 
streighteners; the outer volume computed to 29.95 in. 3 . 

The weight and volume of the mass flow sensor units for the effort in 
subsequent Phases was estimated as follows; Assuming the worst case that 
none of the sensors used in this demonstration phase can be eliminated 
^(realistically, the sensing of pressure and specific heat may no longer 
Jbe needed), the printed circuit boards may be reduced to about 75 x 75 
W in size by a more efficient layout and by using two-sided circuit 
boards. With dimensions in mm, the estimated contributions of the 
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housing and flow channel to volume and weight then are: 83x83x40- 
75 x 75 x 30 + (38 2 - 22 2 )7C/4 x 80, which is equivalent to 22.35 in. 3 
outer volume and 0.995 lb, respectively. The presently used circuits 
weigh less than 0.25 lb. 

The measured power consumption of compensated flow sensors built by 
Honeywell previously was 1.3 watts, when both flow and property sensors 
are operated continuously. The measured response time of the flow sensor 
is about 3 ms, corresponding to a T90 value of 10 ms. A 10-fold flow 
noise reduction benefit may be taken advantage of by letting the T90 
value increase to 1 s. That would leave up to two seconds to perform any 
linearization, temperature and/or pressure corrections; the composition 
correction may not need to be used but once after start-up, or very 
infrequently, depending on the amount of the required correction. 

With over 10 years' experience in the use of such temperature- 
compensated flow microsensors and with over one million of these sensors 
installed, we can look back on a remarkably successful record of 
reliability . In addition, 

1) accelerated life tests indicated that drift rates extrapolated to 
normal operating temperatures and to 10 years of service life would 
amount to less than 0.1% (1) , 

2) direct exposure of the sensor in the proposed flow channel to 
particulate mass loadings of 30 g up to sizes of 300 |im in diameter 
and at mass flows of up to 15 lb/h did not damage the sensor nor 
change its calibration significantly; and 

3) operation and testing of one fully compensated gas meter for 10 
months only altered its flow signal by less than ±1% of reading . 

The above assumes that thermal microconvection effects between the 
sensor elements are negligible, as determined earlier under near- 
atmospheric pressure conditions < 12 > . However, as the pressure increases,, 
the ratio of buoyant/inertial forces, Rbi, increase as well and may 
become large enough to have the associated, thermally-driven flows 
interfere with flow measurement, especially at low flows. The ratio Rbi = 
(AT/T)gL/v 2 = (AT/T) gL (p/p 0 ) 2 /v 0 2 , * shows that this interference is 
proportional to AT, p 2 and L, and inversely proportional to v 0 2 . 


★ 


It may be useful to recap some relevant force groups and their dimensionless 
ratios, to envision the origins of Rbi; see Nomenclature for symbols, except as 


noted: 

- Gravity force, 

- Buoyancy force. 
Inertial force, 

- Surface tension f., 

- Viscous force. 

From these follow, 

- The Reynolds No, 

- The Froude No . , 

The Weber No., 

- The Grashof No., 


GF - 
BF - 
IF - 
SF - 
VF - 
for example: 
, Re - 

Fr - 
We - 
Gr - 


mg, with m « mass 

mgAT/T 

mv 2 /L, 

mcr/ (L 2 p) and 
mvv/L 2 - m(T]/p)v/L 2 . 

Lpv/tl, (-IF/VF), 
v 2 / (gL) , (-IF/GF) or 
v 2 pL/<T, (-IF/SF) . 
L 3 p 2 grT 2 T~ 3 AT, (=Re BF/VF) 


is a hybrid. 


3.3 Actual Design and Fabrication: At the kick-off meeting (8) or 
shortly thereafter it was agreed that: 
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We would use the existing 3x4” PCBs for the demo units, in spite of 
exceeding the dimensional limitations of the hardware we plan to 
replace. This concession may not apply for future production units. 

The output for all gases would be one and the same logio (lb/h) , 
whereby the specified 0-5 V output will correspond to the mass flow 
range from 0.01 to 150 lb/h, for a 15,000:1 dynamic range. 

- The sensor output would provide a digital resolution for NASA's 9-bit 
(2 9 =512) A/D converter of lS^OO 1 / 500 = 1.0194, which is equivalent to 
a resolution of 1.94%, which is still lower than the sensor p-p error 
spec of 6 %, listed as ±3 % in the performance goals of Table 1. 

- We would reduce our Venturi throat area by a factor of 2, so that we 
can meet the external size of the presently used flow channel. This 
will shift the data of Fig.4< 10 > by a factor of 2x to the left. 

- We would achieve this by combining the flow microsensor's differential 

output from 0.01 to 8 lb/h, and its heater-onlv power output from 4 to 
150 lb/h, see Fig. 4 or 7, which show data obtained with a flow channel 
with Venturi throats of 14 and 10 mm ID, respectively. \ 

To meet the performance goals with the 10 mm ID flow channel, we reduced 
the size of the flow sensor header from a T05 to a T018, or from a mesa 
diameter of 0.300" to 0.168", to make it compatible with the reduced 
nozzle diameter, see Figs. 12-14. Under an independently funded 
effort < 14 ), the parasitic thermal microflows were reduced by a factor of 
^870 by a combination of sensor chip design and operation (14) , and 
improvements in the flow channel, which consisted of addition of a flow 
straightener and reduction in the chip head space. Table 2 lists the 
^ised construction materials. 

The electronics, see Figs. 15 and 15a consisted of an analog front end, 
which controlled the sensor heaters to a temperature rise of about 30°C, 
and amplified the sensor element outputs of the flow, temperature and 
thermal conductivity sensor (heater voltage and current) . These were 
then digitized. The thermal conductivity sensor identified the gas being" 
used, so that the stored composition correction could be made without 
requiring any of the usual pressure or specific heat measurements W . The 
communication protocol set up to interface with Carleton is shown in 
Fig. 15b and is self explanatory. 

A photograph of one of the two sensors shipped to Carleton in August of 
1995 is shown in Fig. 16. 

Test results of these units, while meeting expectations in terms of the 
response time, power consumption, insensitivity to orientation up to 200 
psi, and accuracy at flows below 5 lb/h, were disappointing as far as 
accuracies achieved with the high flow sensing approach (above 5 lb/h) , 
labour-intensive adjustments of the matching of the high and low flow 
signals, flow signal errors resulting from the temperature difference 
between the flowing gas and the sensor block, in which the property and 
temperature sensors are mounted, the orientation sensitivity at 900 psig 
and signal noise, and the fact that the remaining orientation 
sensitivity also prevented us from differentiating between N 2 and O 2 . 

these problems were resolved, or recommended that they be resolved as 
follows: 

(1) O 2 and N 2 Differentiation: Output a logarithmic value of volumetric 
flow, so that the differentiation between O 2 and N 2 mass flows (which 
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( 2 ) 


(3) 


(4) 


carries a potential error penalty of 14.2%) is eliminated, except for 
a residual error of less than 0.4%, resulting from the use of 
identical values of Cv in eq. (1) for both O 2 and N 2 . 

S/N Ratio and Output Stability: Upgrade the electronics by taking 
advantage of as many measurements and as much signal averaging as 
possible, which still meets the 3 sec response time; and by 
introducing chopper stabilized amplifiers at the front end, to 
minimize electronic offset drift as possible, 

S/N Ratio, Signal Predictability and Orientation Sensitivity: Provide 
additional averaging by switching to a bypass flow sensing approach, 
see Fig. 17, to enable multiple flow stream pressure taps; the 
implemented design shown in Fig. 17 uses 8 sampling holes for each of 
the two Ap taps and the one property sensor tap. Having the flow 
sensor chip mounted in a well defined microchannel allowed further 
reductions in the orientation sensitivity. In addition, there was 
strong evidence that the heater power sensing mode used at high flows 
could be eliminated, because of the greater range of the more \ 
accurate, differential mode in the bypass design. 

Reduce the thermal mass to the point where the sensor block follows 
changes in gas stream temperature with a shorter time constant; 
thermally isolate both sensors from the sensor housing block; and/or 
incorporate an additional sensor into the flow sensor to enable 
implementation of a more effective temperature compensation. 


Under the same independently funded effort mentioned above < 14 >, such a 
flow sensor was built and made available to this project for test 
purposes. It featured all the improvements listed above except for 
item (4) . 

4. SAFETY AND QUALITY ASSURANCE 



By way of documenting this sensor's fire safety aspects, derived in- 
house before the recent approval for use of this sensor in combustible 
gas environments , we attach ref . (4) in the Appendix. 


The accuracy of our flow calibrator pistons was determined to be within 
a 0.46% standard error margin. This determination included uncertainties 
introduced by dimensions of the glass cylinders, distance between 
optical switches, absolute pressure and temperature measurements, and 
uncertainty in the molar volumes. 


5. LAB TEST RESULTS 


After fabrication and before testing, the sensors were calibrated in a 
process that consisted of a three-step process: 

Determined first were the thermal coefficients of resistance of the 
sensor elements to be used, and the coefficients of the thermal 
conductivity sensor, by operating it over a three-dimensional matrix of 
gases (N 2 , He (or H 2 ) and Ar) , pressures (1-30 bar, limited by our 
pressure sensor) and temperatures (-4 to 60°C) ; then the coefficients 
for the flow sensor were determined in a preliminary fashion at HTC by 
measuring volumetric flows with piston calibrators up to about 15 lb/h, 
at temperatures from about 0 to 50°C, and at pressures up to 900 psi, 
before conditioning the gas down to the ambient levels required by the 
calibrators; these were loaded into the corresponding microcomputer and 
verified. Finally, the units were calibrated throughout the whole range 



/ 
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up to 150 lb/h at Carleton, and an updated set of coefficients was 
■ietermined, reloaded into the microcomputer and verified. 

Fig. 18a illustrates such a verification of Unit#2 calibration for 
thermal conductivity, k. The standard and maximum errors were ±0.28 and 
±0.97%, respectively, with a relative sensitivity of k to errors in 
temperature and heater power of only 1.62, i.e. a worst combination of 
±1% errors would only lead to a ±1.62% error in k. A verification with 
H 2 , He and N 2 is shown in Fig. 18. In both figures, the curves represent 
the true values of k for 3 and 35 bar (45 and 500 psia) . As shown, 
curves for O 2 and N 2 overlap, as shown more clearly in Fig. 19. The 
consequences of this and suggested solutions were already discussed in 
the previous Section. As part of additional sensor checkout, we 
quantified the change in sensor output after a cold start, see Figs. 20 
and 20a, for Units#l and #2, respectively. As shown, the electronics 
takes about 20 min to stabilize; after that time, the readout values of 
temperature-corrected k values has stabilized, as shown, in spite of 
continued increases in temperature. The k values outputs increased by 2 
and 1%, respectively, during warm-up. 

Figs. 7 and 8 presented plots of data to verify the preliminary flow 
calibrations performed before the efforts to eliminate the orientation 
effect was undertaken < 6 * * * * * * * 14 > . Fig. 21 shows the degree of success of those 
efforts, indicating that the orientation effect has been eliminated to 
within pressures of 500 psig, but is still noticeable at 900 psig. 

kThe data in the error-plot of Fig. 22 correspond to Unitfl, after final 
■coefficients were derived and implemented, for the indicated gases, 
pressures and temperatures. As shown, the data fall within the dotted, 
maximum allowed error lines at low flows, but not all do at higher 
flows. Large size points were obtained by the heater-power method, while 
the small-sized points represent those obtained by the differential 
method. The data shown in the lower part of Fig. 22 represent error 
values (right-hand y-axis) of thermal conductivity, Ak in %. 

Flows of He measured with Unit#2 exhibited some unexpectedly high output 
readings, by about a factor of two, even after replacing the flow 
straightener in the flow channel, see Figs. 23 and 24. This cast some 
doubt on our ability to predict its performance with H 2 . Such errant 
behaviour had not been observed with the same sensors when operated with 
microsensors, factory assembled into microchannels, see Fig. 10, or with 
the bypass-based units, see Figs. 25 and 26 (14) . 


6 . CONCLUSIONS AND RECOMMENDATIONS 

The above results indicate that not all performance attributes 

materialized, which had been predicted for operation at pressures up to 

61 bar (900 psi) , based on previous data of flow microsensor performance 

up to 4 bar (60 psia) . The unforeseen factor that increasingly 

interfered as pressure was increased, was one caused by thermal 

microconvection currents generated by the sensor's microheaters. These 

caused both large, signal shifts as well as deterioration of the S/N 
kratio. Whereas the direct influence of these currents could be reduced 
"significantly via implemented improvements < 14 > , the ability of a single 
point flow sensor to meet the accuracy specifications over the large 

15,000:1 dynamic flow range remained marginal. 
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Test results of these improved units, while meeting expectations in 
terms of the response time, power consumption, insensitivity to 
orientation up to 500 psi, and accuracy at flows below 5 lb/h, were 
disappointing as far as accuracies achieved with the high flow sensing 
approach (above 5 lb/h) , labour-intensive adjustments of the matching of 
the high and low flow signals, flow signal errors resulting from the 
temperature difference between the flowing gas and the sensor block, in 
which the property and temperature sensors are mounted, the orientation 
sensitivity at 900 psig and signal noise, and the fact that the 
remaining orientation sensitivity also prevented us from differentiating 
between N 2 and O 2 . 

We therefore recommended that the direct, single-point flow measurement 
of Units#l and 2 be replaced by a bypass, multiple-point flow sensing 
approach. Such recommendation was made on the basis of tests made with 
two such units, which not only showed predictable gas composition N 
correction behaviour, a possibility to replace the heater power approach 
used at high flows by the more accurate differential approach over the 
the full flow range, but also a more linear and predictable output, 
enhanced S/N and the possibility of further reducing the effect of 
orientation-dependent, thermal microconvection, due to the beneficial 
influence of mounting the sensor chip in a well-defined, more 
constraining microchannel . 
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8 . NOMENCLATURE 


c p Specific heat in cal/mol 

Cm Composition correction factor for mass flow, dimensionless 

Cv Composition correction factor for vol.flow, dimensionless 

d Flow channel Venturi throat diameter in mm 

D t Thermal diffusivity in cm 2 /s 

g Acceleration of gravity, in cm/s 2 

G Gauge/sensor output signal, typically in V or mV 
k Thermal conductivity in peal/ (cmKs) 

H Nulled heater power of flow sensor in W 

He Helium gas, gas #2 

Hp Heater power of flow sensor in mW 

H 2 Hydrogen gas, gas #3 

J Nulled flow sensor signal in V 

L Characteristic dimension, in cm or in. 

m Mass flow rate in g/s or lb/h 

M Molecular weight in g/mol 

N 2 Nitrogen gas, gas #1 

O 2 Oxygen gas, gas #4 

p Pressure in bar (=10 5 Pa) or psia 

Pr Prandtl number, Pr = v/Dt = "HCp/ (kM) , dimensionless 
Re Reynolds No., Re = Lvp/rj, dimensionless 
T Temperature in °C 

AT Temperature difference or rise, in °C 
Tr Resistance of temperature sensor 

Ui Heater potential of property sensor in V 

U 2 Bridge potential of property sensor in volts 

v Velocity in cm/s 

V Volumetric flow rate in L/h or ft 3 /h 

rj Dynamic viscosity in microPoises, i.e. pg/ (s cm) 
v Kinematic viscosity in cm 2 /s 

p Density in Stokes, i.e. g/cm3 

a Surface tension 

Subscripts 
c corrected 

d downstream 

H related to heater power 

R related to up- downstream resistance difference 

u uncorrected or upstream 

PBCarl95.W04 ’ 14Aug/NASA 
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TABLE 1. Flow Sensor Requirement Goals*. Rev.: 3 Aug. 1995 


Parameter 

Operating Fluids 

Flow Range, All Gases (lb/h) 

Pressure (psia) 

Max Pressure Drop (psid) 


Requirement 

Nitrogen/Oxygen/Hydrogen/ (He) 
0.01 to 150 
100 to 1000 

0 2 - 0.5 0 900 psig & 15 lb/h 

2.0 0 900 psig & 75 lb/h 

N 2 - 0.5 0 200 psig & 15 lb/h 

2.0 0 200 psig & 75 lb/h 

28.0 0 200 psig & 150 lb/h 


He or H 2 - 0.5 0 200 psig & 15 lb/h 

» 

Gas Temperature (°F) 

Ambient Temperature (°F) 

Ambient Pressure (psia) 


35 to 175 (1.7 to 80°C) N 
40 to 90 (4.4 to 31°C) 

0 to 15 


Acceleration, 

EMI 


vert ical/down 
lateral (g) 


(g) 


£4.4** 

£3 . 3±1 . 5 or £4.8** 
** 


Accuracy, including linearity, 
thermal, pressure and compo- 
sition compensation, and long 

term stability, (%) ±3 of actual flow 0 £0.33 lb/h 


Zero Flow Offset or Error (lb/h) £±0.01 0 £0.33 lb/h 
Volume (in. 3 ) 30** 


Weight (lb) 

Power (W 0 28 VDC) 

Output (VDC Linear Logarithmic) 

Response Time, T90 (s) 

Particle Impact Sensitivity 
Self-Check Feature 
Service Life (years) 


1.35** 

1 . 5 max N 
0 to 5 

<3 

Minimal 

Very desirable 

10, w/o maintenance or recal. 


* Reflecting NASA-JSC's RFO# 9-BE2-11-4-45 req'ts.of 29 Jul.'94 
and those of its Amendment dated 1 Sep.' 94 
** Demo Unit not to be tested on this item; 

N Power will be provided via a 110VAC transformer; the sensor 
electronic components will be rated to 70°C (good commercial 
grade) ; it will implement flow corrections based either on sense 
gas type on-line or receive gas type input via the RS232 . 
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TABLE 2. FLOW SENSOR CONSTRUCTION MATERIALS EXPOSED 
TO FLUID MEDIA. Rev. 14 Aug.1995 


MA T E R IALS* 

Aluminum 6061, anodized 

Aluminum honeycomb, .001" 
Stainless steel mesh, 30-50/in. 


FUNCTION 

Housing 

Flow streighteners 


Glass 

Gold plated on'Ni, on Kovar 
Gold wires, 0.001" diameter 
Rubber (Viton), 0.01-0.03" diameter 
Silicon 

Silicon nitride 
Silicon nitride 


Insulation of header posts 
TO-5 and TO-18 headers 
Connection to sensor chip 
O-ring seals for the TO headers 
Sensor chip 

0.5 to 1 pm thick passivation layer over Pt films 
Sensor heater support membrane, 1 pm thick 


Superglue (Loctitte Superbonder 495) Attachment of Si die to header 




• All at gas temperature, except the <20 mW sensor heater elements of passivated Pt, heated 
to <30°C above gas temperature, which was specified to be <80°C. Max. power only during 
max. flow of <150 Ib/h of N 2 . 
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RG. 1 . PRESENTLY USED ATMOSPHERIC REVITALIZATION PRESSURE CONTROL 
SYSTEM (ARPCS) GAS MASS FLOW SENSOR 
Reproduced from Attm't in NASA-JSC RFO# 9-BE2-1 1-4-45, 29 Jul.'94 
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ICROSENSOR OUTPUTS 
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FIG. 2. MICROSENSOR OUTPUTS AND FLOW CHANNEL 
PRESSURE DROP VS. AIR MASS FLOW. 
NOZZLE THROAT, d=1.4 mm ID. 

WINDTUNNEL AIR AT AMBIENT CONDITIONS. 
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d, THROAT DIAMETER IN mm 


1. FLOW VELOCITIES IN THE VENTURI THROAT 
AT 75 Ib/h OF N 9 , VS. THROAT DIAM. 
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FIG. 2a. RESOLUTION OR GRAININESS OF A 9-BIT A/D PROCESSOR VS. MASS FLOW, 
DEPENDING ON LINEAR OR LOG SENSOR SIGNAL OUTPUT. 

THE CHOSEN LOG OUTPUT FOR THE RANGE FROM 0.01 TO 150 Ib/h PROVIDES 
A FLOW-INDEPENDENT RESOLUTION OF ±1.94%. 
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FIG. 1. RESOLUTION OF A DIGITAL PROCESSOR AS 
A FUNCTION OF ITS WORD SIZE IN BITS 
AND THE DESIRED DYNAMIC RANGE 
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FIG. 9. FLOW MICROSENSOR DIFFERENTIAL OUTPUT VS. FLOW AT VARIOUS 

PRESSURES, SHOWING AN INCREASING NATURAL CONVEQTION EFFECT WITH 

DDCCOt IDC 
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FIG. 5d. MBRIDGE SENSOR OUTPUT CORRECTED FLOW 
ERROR FOR SENSOR NO.: #2347, USING 
CORRECTION DERIVED FOR SENSOR #2348 



























SENSOR CIRCUIT: ANALOG FRONT-END FOR k, T and V 
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Flow 

Sensor 


Property 

Sensor 


Hp » V UiH (Hp)* » m C H 

AR(j U =Gdu » Vu ( R (ARdu) » mcR 

null 


100 ms 900 ms 100 ms Cm 


m » x=(5/4.176091)logio(100m) 
(5/9.61 5806)log e | 

I 

* 

Output Pin #6 
0<Sxin log(lb/h) £ 5 V 


Ui,U 2 I 

null, TR»T,k,gas- J 


100 ms 900 ms 100 ms 


I/O is connected via the following pin numbers of the 13-pin connector: 
-INPUTS OUTPUTS 


1. +5 V Digital 


2. Digital Ground 


3. +1 5 V Analog 


4. Analog Ground 


5. -15 Analog 

02, N2, H2, He, 03 or N3 

6. 

D/A Output -x / 

7. 

D/A Return f 

8. 

RS232 Transmit - U#,X,T,Ui,U 2 , gas.Gdu.Hp.m CR LF 

9. 

RS232 Receive - 02, .... or N3 

10. 

RS323 Common 

11. 

Digital Ground 

12. 

Analog Ground 

13. 

Chassis Ground 


Rev.: 21 Auq.*95 


* The corresponding coefficients are upgradable on Serial EE PROM, 24C32 
** Data string transrnitted once/second; ASCII code: ±x.xxxxxxxE±xx for each of the listed 
variables, which are comma-separated and of which each set is finished with a carriage 
return (CR) and a line feed (LF); received inputs: any of the first four X above for 

individual gases, but enter "N3" for N 2 /He/H 2 , or "03" for 02/He/H 2 ; the verification 
-response is the same X as part of the transmitted string. 

FIG. 15b. MASS FLOW SENSOR #2: SIMPLIFIED MEASUREMENT PROTOCOLAND 
INFORMATION FLOW CHART, LEADING TO ANALOG 0-5 VDC OUTPUT. 

Both sensors are pulsed at about 100 ms; the gas is reidentified during each cycle; 
the corresponding flow correction factor, Cm, is recomputed and reapplied to the flow 
sensor signal. The flow signal is downselected from two measurements as shown. 
The phosen mass flow, m, is then converted to its logarithm and converted to analog. 




FIG. 16. PHOTOGRAPH OF SMART MASS FLOW SENSOR DEMONSTRATION UNIT, 
SHOWING THE IN- AND OUTLET FITTINGS IN THE LOWER LEFT BLOCK. 


4.0“ 


Electronics 


LSSSl 


X99005-AW 


U5iJ 


FIF. 4a. IMPROVED FLUID FLOW SENSING SYSTEM BASED ON A BYPASS ACROSS A 
LAMINAR FLOW RESTRICTION (CENTRAL HONEYCOMB). FLUID REACHES 
BOTH PROPERTY AND FLOW SENSORS DIRECTLY. A METAL BRACKET HOLDS 
THE FLOW SENSOR IN PLACE AGAINST FLUID PRESSURE, AND ITS LEADS 
ARE ACCESSIBLE W/O ADDITIONAL FEEDTHROUGHS. 
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FIG. 5. WARMUP TIME OF COMPENSATED FLOW 
SENSOR W/C0NST.AT htr =30°C DRIVE. 
START FROM COLD ELECTR. & SENSOR. 
0#5 FILE: NS224N01 .CSV, p = 1 bar, NO INSUL. 
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FIG. 26. MASS FLOW MEASUREMENT W/MICROSENSOR, 
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FIRS SAFETY ASPECTS OF Si3N 4 -Pt -BASED 
MICR0BRZD6E FLOW SENSORS 

by 

U. Bonne, T. Stratton and R. Horning 
Honeywell Sensor and System Development Center 

May 20, 1992, rev. May 29, 1992 


INTRODUCTION 

The fire safety and flammability of Honeywell's flow micro- 
sensor (1*2,3) described here from two points of view: 

The hazard resulting from operation in a flammable gas such as 
natural gas-air or hydrogen-air mixtures and 

K 

The hazard resulting from operation in a pure, high pressure 
oxygen environment. 

These hazards could be thought of, in the worst case, to lead to 
fires involving the sensor elements as ignitors of either the 
flammable mixture or of the sensor hardware itself. 

To determine the severity of these hazards we checked whether: 

1. The sensor heater, during normal operation, could ignite a 
surrounding gas mixture or 

2 . The sensor heater could ignite itself during normal operation, 
to burn in an oxygen atmosphere or 

3 . The sensor heater could ignite a gas mixture under abnormal 
operation or 

4 . The sensor heater could ignite itself during abnormal operation 
to burn in an oxygen atmosphere. 

5. Any part of the sensor or its flow channel hardware could 
support combustion in air or in an oxygen atmosphere, if 
ignited by an external source. 

The difference between normal and abnormal operation is defined, 
for the purposes of this review, by the difference in heater 
power: under normal conditions the approx. 50 x 100 Hm sized 
heater is raised 60-1 60°C above ambient temperature, with about 2-6 
mw of input power, and by applying less than 10V DC; under 
abnormal conditions, an over-voltage may first overheat and 
eventually vaporize the platinum thin film, and cause an electric 
spark and an open circuit . 


DISCUSSION 


In addressing the above five hazard modes, the physical evidence 
and experience which lead us to believe that this sensor is safe, 
is as follows: 

1 . Auto ignition temperatures of hydrocarbon-air mixtures range 
from 537°C for methane, 466°C for propane, to 206°C for n- 
nonane/n-hexadecane (4) . At these temperatures, large volumes of 
stoichiometric gas mixtures will eventually ignite. However, 
ignition by wires at these temperatures is much more difficult: 
a bare nickel or platinum wire needs to be heated to over 1000 
and 1150°C, respectively to achieve natural gas-air ignition 
under lean conditions, and still higher temperatures under 
stoichiometric or rich conditions. Given the smallness of the 
microbridge heater (i.e. well below even the quenching 
distance) and its normal operating temperature of less than 185 
to 285°C, even under 125°C-under-hood-automotive ambient 
temperatures, ignition is extremely unlikely if the heater is 
operated at a peak absolute temperature of less than 200°C. 

2. Reaction rates for Si 3 N 4 to form Si02 and N 2 or NO 2 are so slow, 
even at 500°C, that combustion is not sustained. Evidence: 

a) The most energetic reaction, forming 3Si02 + 2 N 2 , releases 

473 kcal per mole of Si 3 N 4 < 6 ) (465 kcal to form 3S102+ 2 NO 2 ) , or 
158 kcal/mol of Si02 formed, which is less than released by the 
reaction: Si + O 2 » Si02 + 209 kcal/mol. 

b) The Si-oxidation reaction can also not sustain itself, 

because its heat release rate is several orders of magnitude 
lower than the rate of heat dissipation by thermal conductivity 
into the surrounding oxygen atmosphere: The known "thin oxide 

layer" or reaction-limited oxidation rate of Si, at 500 *C, 
through its adhering oxide, is 3 .5*10" 7 Jlm/h, which is equivalent 
to 6*10 -7 cal/(h cm 2 ). The heat conduction rate, assuming only a 
1000 *C/cm gradient, is over 8 orders of magnitude larger. 

c) The reaction-limited oxidation of Si 3 N 4 at 1200 *C was 
computed to be 0.052|lm/h^ 11 ^ , which still is equivalent to a 
factor of lOOOx less than the heat dissipation rate. 

d) Even if oxidation occurs at high oxygen concentrations and 
temperatures, the forming oxide adheres to Si 3 N 4 and Si02/ as 
with the oxides of aluminum and stainless steel. 

3. One abnormal condition of operation was tested. A microbridge 

flow sensor was introduced into a 0.8 stoichiometric natural 
gas-air mixture. An over-voltage applied to the heater, causing 
instant vaporization and a little spark; but no ignition of the 
combustible mixture resulted < 7 ) . The applied heater voltage was 
lOx the normal level required for slow burnout, but was not 
sufficient to generate enough energy (10~ 5 J) to exceed the 
minimum ignition energy of 2*10“ 4 . 

The auto ignition temperatures of hydrocarbons with oxygen are 
only about 5°C lower than with air, except for butane (about 
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130°C lower) and possibly propane (not listed in ref. 4), but 
none is lower than 200°C. 

4. The reasoning in above item 2 also holds here: while an over- 
voltage may disable the microbridge element (heater or sensor) , 
the achieved local heating will not support a propagating 
combustion wave to reach the Si chip. Solid Si is not known to 
burn in an O 2 environment, unless ground and dispersed as a 
cloud of particles. Then its heat of combustion would be 209 
kcal/mole . 

5 . The most vulnerable part of the exposed part of the sensor 
would be the honeycomb and screen (s) used to condition the 
flow. Therefore their composition was selected to not be 
flammable even in an oxygen atmosphere. Aluminum alloys 11XX- 
ANN (MIL-A-148) or 2024-T81 in foil< 10) and 1100 or 2024-T4/T6 
in wire form are rated for use in oxygen up to lOOOpsi in 
Rockwell's Materials Directory. Other metals such as braze Cu 
and Ni, Ni-Cr-Fe wire, and Weld-wire Ni-based are also listed .. 
with ratings of up to 3000psi in oxygen and are candidate 
materials for screens. 


The used 25 Jim OD (0.001") gold wires to connect the sensor chip 
with the TOS header are not viewed as a hazard in an oxygen or air 
environment. 

Fused zener diode barrier circuits can be added to the sensor 
input if a probability exists for the accidental application of 
voltages above those handled safely by the circuit (typically 
30V), thus redering the circuit "intrinsically safe" (8) . The 
"internal" sensor and circuit are intrinsically safe since its 
maximum heater energy levels (<lmA x <15V) , in case of 
malfunction, would be more than lOOOx below the level needed to 
ignite combustible mixtures of hydrogen-air ( 8 ) , in the presence of 
Al, Cd, Mg or Zn, or even of hydrogen-oxygen mixtures, with about 
lOOx lower ignition energies, if the measurements with 
hydrocarbon-air and oxygen may be used as a guide. 


CONCLUSIONS 

In conclusion, whether operating normally or under extreme 
(accidental) conditions, the hazard of having the Honeywell 
Si 3 N 4 /Pt/Si flow sensor self-ignite or cause a combustible gas-air 
mixture to ignite is so small that it has not been measurable to 
date . 
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i GAS FLOW SENSORS FOR SPACE APPLICATIONS 

f Phase I. Fabrication of Two Demonstration Units 

1 . EXECUTIVE SUMMARY 

The purpose of this program was to demonstrate the performance of a 
thermal flow microsensor. It would be intended for the measurement of 
flows of helium, nitrogen, oxygen and hydrogen in both the ARPCS 
(Atmospheric Revitalization Pressure Control Subsystem) and in the FCS 
(Fuel Cell System) , which have been using sensors that experienced high 
failure rates. 

An existing sensor design was adapted for this project. Two sensors were 
built and tested, and delivered to Carleton for further tests, after 
improvements were made to overcome interferences by what was recognized 
as thermal microconvection. 

\ 

In spite of* the improvements, the test results indicate that not all 
performance attributes materialized, which had been predicted for 
operation at pressures up to 61 bar (900 psi) , based on previous data of 
flow microsensor performance up to 4 bar (60 psia) . The unforeseen 
factor that increasingly interfered as pressure was increased, was one 
caused by thermal microconvection currents generated by the sensor's 
microheaters. These caused both large, signal shifts as well as 
deterioration of the S/N ratio. Whereas the direct influence of these 
currents could be reduced significantly via implemented improvements , 
the ability of a single point flow sensor to meet the accuracy 
^specifications over the large 15,000:1 dynamic flow range remained 
marginal . 

Test results of these improved units, while meeting expectations in 
terms of the response time, power consumption, insensitivity to 
orientation up to 500 psi, and accuracy at flows below 5 lb/h, were 
disappointing as far as accuracies achieved with the high flow sensing 
approach (above 5 lb/h) , labour-intensive adjustments of the matching of 
the high and low flow signals, flow signal errors resulting from the 
temperature difference between the flowing gas and the sensor block, in 
which the property and temperature sensors are mounted, the orientation 
sensitivity at 900 psig and signal noise, and the fact that the 
remaining orientation sensitivity also prevented us from differentiating 
between N 2 and O 2 . 

We therefore recommended that the direct, single-point flow measurement 
of Units#l and 2 be replaced by a bypass, multiple-point flow sensing 
approach. Such recommendation was made on the basis of tests made with 
two such units, which not only showed predictable gas composition 
correction behaviour, a possibility to replace the heater power approach 
used at high flows by the more accurate differential approach over the 
the full flow range, but also a more linear and predictable output, 
enhanced S/N and the possibility of further reducing the effect of 
orientation-dependent, thermal microconvection, due to the beneficial 
influence of mounting the sensor chip in a well-defined, more 
^constraining microchannel . 
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Figure 4-31. Curve of resistance circuits, cadmium disc. 
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